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Abstract. The failure of materials and interfaces is mediated by cracks, nearly 
singular dissipative structures that propagate at velocities approaching the speed of 
sound. Crack initiation and subsequent propagation - the dynamic process of fracture 
- couples a wide range of time and length scales. Crack dynamics challenge our 
understanding of the fundamental physics processes that take place in the extreme 
conditions within the nearly singular region where material failure occurs. Here, we 
first briefly review the classic approach to dynamic fracture, "Linear Elastic Fracture 
Mechanics" (LEFM), and discuss its successes and limitations. We show how, on 
the one hand, recent experiments performed on straight cracks propagating in soft 
brittle materials have quantitatively confirmed the predictions of this theory to an 
unprecedented degree. On the other hand, these experiments show how LEFM 
breaks down as the singular region at the tip of a crack is approached. This 
breakdown naturally leads to a new theoretical framework coined "Weakly Nonlinear 
Fracture Mechanics" , where weak elastic nonlinearities are incorporated. The stronger 
singularity predicted by this theory gives rise to a new and intrinsic length scale, 
t n \. These predictions are verified in detail through direct measurements. We then 
theoretically and experimentally review how the emergence of l n \ is linked to a new 
equation for crack motion, which predicts the existence of a high-speed oscillatory 
crack instability whose wave-length is determined by i n \. We conclude by delineating 
outstanding challenges in the field. 



PACS numbers: 46.50.+a, 62.20.mm, 62.20.mt, 89.75.Kd 
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1. Introduction 

Our fundamental understanding of why and how materials break is, surprisingly, only 
about 100 years old. Estimates of the theoretical strength of materials range from 
E/tt to E/8 PQ, where E is the Young's modulus of the material. These estimates 
are simply related to the work that one would have to expend in order to overcome 
the potential well that holds neighboring atoms together. When pulling on a sheet of 
window glass, for example, one would therefore expect it to stretch at least 10% before 
breaking. In practice, if you are very careful, you might be able to reach a hundredth 
of that strain before the glass breaks. This huge disparity between the theoretical and 
practical strengths of materials is entirely general in brittle materials. Obviously there 
is something that these calculations are missing. The answer to this is the existence of 
cracks in these materials. 

The modern history of fracture mechanics started with a calculation by Inglis in 
1913 [2j. Inglis found that the imposition of an elliptical hole into a linear elastic sheet 
under uniformly applied tensile stress, entirely changed how stresses are distributed in 
the vicinity of the hole; externally applied stresses are amplified at the tip of the larger 
axis by the ratio of the large and small axes of the ellipse. This stress amplification 
increases without bound if the ellipse is "squashed" to form a crack. In this limit, 
the stress tensor at a crack's tip, <r, becomes singular, increasing as a ~ 1/v^ 7 where 
r is the distance from the tip. This singular behavior at a crack's tip is the basis 
for understanding both the strength of materials and the physics of fracture. In the 
example of the glass plate a crack of length 10/xm is sufficient to reduce the plate's 
theoretical strength by a factor of 100. Simply stated, material strength is governed by 
the formation and subsequent propagation of cracks. A material fails when the tip of a 
crack at its weakest point starts to propagate. 

Models describing the propagation of cracks, or dynamic fracture, only date back to 
the second world war. In fact, one of the first theories of crack propagation was derived 
by Sir Neville Mott who was enlisted to understand the causes of brittle fracture [3] . A 
specific problem of interest was why a significant number of rapidly constructed cargo 
ships, known as the Liberty ships, underwent cataclysmic failure either immediately 
upon or shortly after their initial launch. Since then, an immense amount of research into 
how and when cracks propagate has been performed. Despite this effort, there remain 
significant and fundamental aspects of crack propagation that we do not understand. 
We will show, in this review, that an improved understanding of the material behavior 
in the vicinity of the singular region surrounding a crack's tip has shed a new light on 
a number of these aspects. 

1.1. Why is fracture interesting to a physicist? 

Understanding fracture is clearly important in practical applications such as designing 
stronger materials and structures. The process of fracture also raises a variety of 
important and interesting physical questions. Crack propagation involves the integration 
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of essential physics at an extreme range of widely varying spatial and temporal scales 
that link the macroscopic scales where energy is injected, stored and transported to the 
microscopic scales where dissipation in the form of irreversible material deformation and 
fracture takes place. Fracture processes therefore couple a huge range of scales. At the 
smallest scales one might wish to understand how the introduction of nano-structures 
can affect the ultimate strength of a material. At geophysical scales one would like to 
understand what scales are needed to determine whether a natural fault will lose its 
stability and generate a massive earthquake. The study of fracture dynamics has also 
brought to the forefront numerous fundamental questions that pit continuum theories 
against discrete physics. Questions of when and how do atomic scales come into play 
can be important [4J [5j [6] . 

A crack is characterized by the singular stress fields that drive it. As nature 
generally will not allow "real" singular stresses to develop, the toughness of materials 
is essentially determined by how these singularities are regularized. The mechanisms 
by which nature performs this regularization determine whether a given material can 
be used as a structural material (e.g. steel) or not (e.g. window glass). Furthermore, a 
propagating crack can rapidly reach velocities that approach material sound speeds, the 
speeds at which information propagates in these systems. Thus, the physics of crack 
propagation are closely related to questions of the formation and regularization of finite 
time singularities that are formed at moving fronts. These fronts are "relativist ic" in 
the sense that they travel at speeds approaching the information speed in the material. 

Achieving a fundamental understanding of the dynamics of fracture may also 
shed light on a rather broad class of conceptually related physical systems. Fracture 
propagation is a rather close relative of the broad class of physical problems that 
can be loosely characterized as "growth" problems. These systems involve situations 
where a moving boundary separates two distinct phases. Examples include propagating 
fronts as vehicles for phase transitions (e.g. the spread of thermal convection fronts 
or the motion of an interface between stable and unstable regions) [71 [8], Laplacian 
growth problems (e.g. crystal growth, Saff man- Taylor like problems of fluid invasion 
or imbibition) [9j, the physics of lightening [TO] , flame propagation pj] and reaction- 
diffusion fronts, and the general problems of interface propagation and roughening (e.g 
KPZ-like problems) [12]. A common denominator of all of these problems involves 
coupling of two media, whose behavior is described by (often linear) field equations (e.g. 
Laplacian or diffusion equations, wave equations, Ginzburg-Landau equations) that are 
coupled at the boundary between the two media [9j. This boundary (or front) generally 
undergoes rich and varied space-time dynamics that are not a priori known, but are 
rather determined self-consistently from the global solution. Precisely these dynamics 
are what we would like to understand. In a 2D crack propagation problem, the boundary 
is the crack tip itself and the two crack lines left behind it. In a 3D body a crack is a 2D 
sheet whose leading edge is a putatively singular one-dimensional front that separates 
intact and fractured material. The dynamics and instabilities of this rapidly propagating 
singular front result from coupling the space-time behavior of the intact medium, as 
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described by linear wave equations, to the moving boundary, which is defined by stress- 
free boundary conditions at the crack faces. Via the surrounding fields, the front can 
interact in space and time with both itself and the crack surface that it had previously 
formed [T3l [Tlj . In addition, the crack front dynamics can be affected by interaction 
with material inhomogeneities in its path [131 HH H51 HEl HZl UHl HSl [20l EH [22l [23] . 

The study of dynamic fracture is also an interesting example of a physical system 
that is strongly out of equilibrium. There are few better examples for this than a 
nearly singular tip that is pulling a material apart at nearly the speed of information 
propagation. This challenging problem has forced us to extend our knowledge of 
materials well past the point where equilibrium properties are useful. An instructive 
example of this was encountered by Marder and coworkers (4[ [24] when they incorporated 
potentials that successfully describe equilibrium properties of silicon [25j in a molecular 
dynamics simulation of fracture in this material. Using this potential, cracks were barely 
able to propagate and the energy needed to initiate fracture was 400% larger than the 
measured value |H[26]. Later work suggested that, in possibly the most widely studied 
material in history, quantum mechanical calculations were needed to enable an accurate 
quantitative description for the simplest mode of crack propagation (the problem of a 
single propagating crack) [271 EE] • 

Modern theories of amorphous plasticity also have their roots in our collective 
frustration in finding a fundamental understanding of the onset and propagation of 
cracks [29], [30]. Entire classes of continuum models [HU [32l [33], [34] that purport to 
describe the generic breakdown of the singular behavior at a crack's tip in amorphous 
materials have been shown to be unable to describe how and why a simple propagating 
crack becomes unstable. Failed attempts to do so have lead a number of groups 
to initiate theories to provide a better and fundamental description of how plastic 
deformation takes place in amorphous materials, where plasticity due to dislocation 
formation and propagation cannot play a role [29], [30], EU [35j [36], EH [38j [39] . 

1.2. The breakdown of linear elasticity near crack tips 

Naively, one would think that the existence of a singular stress at a crack's tip should 
be sufficient to cause the material at the tip to fracture. Assuming that Hooke's law 
describes the elastic behavior of a material (i.e. that materials are linearly elastic - 
stresses are proportional to strains), singular stresses should lead to singular strains, 
hence to bonds at the tip of a crack being pulled apart as a crack continually extends 
itself. In practice, there is a threshold for crack propagation, which in essence defines 
the toughness of a given material. Where does this threshold come from? 

The answers to this question are related to the ways that nature manages to 
regularize this stress singularity. Obviously, stresses and/or strains cannot really be 
mathematically singular so something has to happen as one approaches the small length 
scales that are near a crack's tip. The region where the 1 / yjr singularity of the stress field 
breaks down has been named the "process zone" . Within the process zone all dissipative 
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and/or nonlinear processes that take place when materials are stretched beyond the 
point where linear elasticity breaks down are assumed to take place. The process zone 
is, essentially, the rug under- which we sweep all of the "dirty" processes for which 
we either lack fundamental understanding or where continuum theories break down. 
This nebulous region encompasses processes such as plastic (irreversible) deformations, 
nonlinear elastic effects, damage accumulation, visco-elastic processes, and dissipative 
effects due to discreteness at the atomic scale. 

What is the size of the process zone? There are few direct measurements of this 
region and its size estimates vary over a wide range from material to material. Rough 
estimates are sometimes obtained by equating the singular field to a material's yield 
stress. These estimates are often very crude and can vary considerably. In some 
materials the natural cutoff for the singularity is the size of the discrete atomic scale 
where the continuum theory that gave rise to the stress singularity must certainly break 
down. In brittle materials such as crystalline silicon and even amorphous glasses nm 
scales have been quoted for this region [32j. In brittle acrylics estimates range from 1- 
10/iTO. In very tough materials such as aluminum the process zone can reach mm scales 
[32] , as the size of the dissipative zone is determined by the scale at which dislocations 
moving out of the crack tip pile up and lose their mobility. 

A key assumption of fracture mechanics that provides a way to circumvent the 
myriad material-dependent dissipative mechanisms and nonlinear processes that are 
purported to take place around a crack's tip is the assumption of "small-scale yielding" . 
Small-scale yielding assumes that the details of the structure of the process zone can 
be largely ignored, if the region where these processes take place is sufficiently small 
compared to other characteristic length scales in a given problem. In such a case it 
is possible to justify a separation of scales; stress fields outside of the process zone 
can be described to high accuracy by the singular contribution to the linear elastic 
stress field, a ~ K/y/r which, at intermediate scales, dominates all non-singular linear 
elastic contributions. This separation of scales enables one of the triumphs of fracture 
mechanics: a universal description of the functional form of the linear elastic stresses and 
strains surrounding the tip of a crack. As energy dissipation is confined to a small region 
within this universal field, predictions of material failure and the consequent motion of 
a crack are entirely determined by the value of X, which is called the "stress intensity 
factor" . Knowledge of K is equivalent to knowing the energy flowing into to the crack's 
tip. The assumption of small-scale yielding suggests an explanation for the universality 
of fracture behavior; e.g. why does a brittle plastic fracture in the same way as the glass 
in your window - despite that fact that the dissipative processes in both materials are 
vastly different. This implied universality, which is sometimes termed "K-dominance" , 
suggests that all one needs to know is how to either measure or compute K for a given 
external application of stresses in order to predict how a crack in a given material will 
behave. 
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1.3. Linear Elastic Fracture Mechanics (LEFM) 

Linear Elastic Fracture Mechanics, or LEFM, provides the basis for our current 
understanding of fracture j32j [33l [40j EJ S21 [43l |44j. In this section we will both 
describe LEFM and set up the notation that will enable us to later describe fracture 
mechanics when the assumptions of linear elasticity are relaxed. 

We start by defining the motion which is assumed to be a continuous, 

differentiate and invertible mapping between a reference (undeformed) configuration 
described by coordinates x and a deformed configuration described by coordinates x\ 
both in K 3 , 

x f = 4>(x, t) = x + u(x, t) , (1) 

where u(x,t) is the displacement field. The components of the displacement gradient 
tensor H are defined as 

Hij = djUi . (2) 

Note that partial spatial derivatives are assumed to be taken with respect to the reference 
configuration x ) unless otherwise stated. However, the distinction between the reference 
and deformed configurations manifests itself only to nonlinear orders in H and hence 
it makes no difference in the framework of linear elasticity. Within this framework one 
assumes \H \ <Cl and defines the infinitesimal strain tensor as 

e^liH + H^, (3) 

where the superscript 4 T" stands for the transpose of a tensor. The linear elastic energy 
density functional U (s) of isotropic materials can be expressed as 

U(e) = -\{tref + litre 2 , (4) 

where A and \i are the Lame constants [45]. The Cauchy stress tensor a is 
thermodynamically work-conjugate to s, cr = d s U(e) = \ (tre) I+2fie (I is the identity 
tensor), which is nothing but Hooke's law [45j. Substituting the latter in the linear 
momentum balance 

V • <r = podttu , (5) 
where p is the reference mass density, we obtain the standard Lame equation 

/iVV 1 ) + (A + ^)V(V • u^) = Po d tt u^ , (6) 

which is the basic equation of linear elasticity. The superscript (1) was introduced to 
stress the fact that this equation is a first (linear) order approximation in the magnitude 
of H. This will be important later when higher order contributions will be discussed. 
Note that angular momentum balance is automatically satisfied due to the symmetry 
of the Cauchy stress, cr = cr T [45J. 

Consider now a long straight crack propagating in a 2D body and define a fixed 
Cartesian coordinate system (x, y) such that the crack's tip propagates steadily at a 
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velocity v in the positive x-direction. The ^-direction is perpendicular to the crack's 
faces. A crack is defined physically as composed of surfaces that cannot support stresses, 
i.e. by the following boundary conditions 

(J xy (r,Lp = ±it) = (j yy (r,(p = ±ir) = . (7) 

(r, (p) is a polar coordinate system that moves with the crack tip, which is related to the 
rest frame by r— yjjx — vt) 2 + y 2 and (^ = tan _1 [?//(x — vt)]. Mathematically speaking, 
a crack can be regarded as a moving branch cut. Under steady-state propagation 
conditions we expect all of the fields to depend on x and t through the combination 
x— vt and therefore d t — —vd x . tp = coincides with the positive x-direction and Lp = ±n 
define the two opposite crack's faces. Using Hooke's law, equations ([7]) can be rewritten 
as 

— jir^d^u^ — jidrdjp = , 

- (A + 2/x)r- 1 d^ 1) - A d r u^ = , (8) 

for if ±7T. 

For mode I (tensile) symmetry, i.e. u^t\x,—y) = u x l \x,y) and u^\x,—y) = 
— in either plane-strain or plane-stress conditions j33j[43], the two-term near 
tip asymptotic solution of equation ^ with the boundary conditions of equations ([8]), 
is given as [33j 

,(D = ™* 



u y 



(i) _ _ 2Ki®d 



(! + a l) r d 2 cos (y ) - 2o^ s r s 1/2 cos (y 



(A + 2fi) T r cos 



/xV27r£>0) 



H ere a 2 d s = 1 - v 2 /c\ s tan y? d>s = a d , s tan (p, r d;S = ry/l - (v sin (f/c d ^ s ) 2 . c d = 
y/(X + 2/x)/po and c s = ^//x/po are the dilatational and shear wave speeds, respectively. 
Finally, Kj is the mode I "stress intensity factor" , T is the "T-stress" (the amplitude of 
the sub-leading term in the linear elastic asymptotic expansion) and D{v) = Aa s a d — (1 + 
a^) 2 . The latter vanishes at the Rayleigh wave-speed cr, D(v = cr) = 0. Note that the 
term proportional to Kj in equations Q gives rise to the famous 1 / ^Jr displacement- 
gradients (and stress) singularity discussed above [33j [43] . Similar expressions are 
obtained under global shear loading (mode II fracture), where the mode II stress 
intensity factor Kjj appears instead of Kj. The role of Kjj in determining the direction 
of crack propagation will be discussed later. 

One immediate implication of the solution in equations ^ is that the crack that 
is represented by a straight branch cut in the reference configuration becomes parabolic 
near its tip in the deformed configuration 

<P x (r,±n) = -Xi<Pl(r,±n) , (10) 

obtained by expressing x (r, ±7r) = — r + u x (r, ±7r) as a function of y (r, ±7r) = u y (r, ±7r) 
(see equation ([T]) for the definition of 0), where the curvature xi can be easily read off 
equations ^ in terms of Kj and T. This is the so-called parabolic "crack tip opening 
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displacement" (CTOD). It is important to note that the values of Kj and T in equations 
(|9j)-( 10 ) depend on the driving of the system, specifically the globally applied stresses. 



These, of course, cannot be determined by asymptotic analysis. In the next subsection 
we will see how the stress intensity factor is related to both the applied stresses and the 
energy flux into the crack tip, and consequently to the crack's motion. 

1.4- Fracture initiation and energy balance 

When a crack extends, it forms two new surfaces at an energy cost. We define the 
fracture energy, T, as the energy needed per unit extension of a crack to create these 
surfaces. Griffith suggested "energy balance" as a criterion for crack stability [46]. The 
resulting "Griffith criterion" states that a crack will lose stability when the change in 
potential energy U in the surrounding medium (including the loading machine) released 
by the crack upon an infinitesimal extension, 5/, of its length / surpasses T. While 
the "Griffith condition", dU/dl>T, avoids the need to explicitly take into account the 
nature of the stresses at a crack's tip, it is clear that the energy needed to propagate a 
crack and the form of the stress field surrounding a propagating crack must be related. 
(We note in passing that Griffith considered only the bare surface energy 2j as the cost 
of crack initiation. The fracture energy T is a generalization of the concept of surface 
energy, see below for an additional discussion of this point). As we have seen in equation 



(10), sufficiently close to the crack's tip, the only part of the asymptotic solution that 
can contain information about the global loading conditions is the stress intensity factor, 
Kj. Considering a nearly static crack, Irwin indeed showed that the energy flowing into 
the crack's tip per unit crack extension, G(v), is given (under plain-strain conditions 

m) by m 

GUo = l ~^K 2 i • (11) 

In the presence of shear stresses near the crack tip an additional contribution 
proportional to K]j appears [33j . 

The quantity G in equation ( pTj ) is called the energy release rate (even though it 
involves no rate, its dimensions are energy per unit area). It quantifies the amount 
of energy per unit fracture surface that is flowing into the tip of a crack, where it is 
dissipated. The generalized Griffith condition for fracture initiation therefore reads 

1 - v 2 2 

G v ^o = — ^ — Kj = T v ^ . (12) 



Equation (12) demonstrates that the generalized Griffith condition for the stability of a 



static crack is equivalent to positing a critical value of the stress intensity factor. The 



dimensions of Kj, as can be seen in equation (|9j), are those of stress x ^/length. The 
explicit value of Kj can be calculated for a static crack in a given loading configuration. 
For example, when a constant tensile stress is applied to the remote boundaries of 
a large sheet in which a crack of length I exists, Kjoca^ VI. If the same uniform stress 
is applied to an infinitely long strip of width 2b then Kj oc Vb. 
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Whereas both Kj and G are loading-dependent quantities, T is considered to be 
a material-dependent quantity that simply specifies a given material's resistance to 
being broken, Y(v) can in fact be a rate-dependent function and is, therefore, generally 
dependent on the instantaneous crack velocity v. Y{v) encompasses all of the dissipative 
processes that take place within the process zone. It is not simply the energy cost of 
breaking a single plane of material bonds, i.e. the surface energy 27, as a variety of 
possible dissipative material-dependent processes need to take place to enable a bond 
to be finally broken. Often these processes require orders of magnitude more energy 
that the simple fracture of the bonds that hold a material together. For example, 
the fracture energy of a brittle acrylic such as PMMA (poly-methyl-methacrylate or 
"plexiglas") requires 500 — 1000 J/m 2 whereas breaking a plane of dense carbon bonds 
would typically require about two orders of magnitude less energy. The huge amount of 
"extra" energy that goes into Y for PMMA is thought to result from plastic deformation 
of the polymer that is a necessary condition for the separation of the polymer chains 
that compose PMMA. 



1.5. Equations for crack growth rate 



Equation (12) tells us what is needed to initiate fracture. What happens once a crack 
starts to propagate? The condition of dynamic energy balance, G{v) = Y(v), which 
further extends the generalized Griffith criterion to all velocities (crack growth rates 
v) is the basis for the formulation of a dynamic theory of fracture [33] . Once G(v) 
is calculated and Y(v) is either calculated or measured, one is able to predict v as a 
function of variables such as the crack length and parameters such as any externally 
applied stresses as long as the crack path is known a priori. That means that the scalar 
equation G(v) = Y(v) is, in principle, capable of determining the crack growth rate but 
not the crack's direction of propagation. The analytic calculation of G(v) is generally 
very difficult. Two crack configurations for which it has been performed are: 

• A medium of effectively infinite spatial extent for any combination of forces applied 
to the crack faces (or any problem that can be mapped to such a problem) [33j . 

• A semi-infinite crack propagating at any velocity, v, within an infinitely long strip 
whose finite boundaries are subjected to constant displacement conditions [47J. 

Both of the above calculations are limited to simple, perfectly straight cracks that 
have not undergone any path instabilities, such as oscillations or branching. Propagating 
cracks in each of the above calculations have the same form of asymptotic stress tensor 
field, a oc Ki{v)/y/r, but with different dependencies of Kj(v) in terms of the system 
parameters and geometry. Below we will briefly review the main features of each of the 
resulting equations of crack growth rate. 



1.5.1. The motion of a crack in an infinite medium The first calculation is due to 
Eshelby [48], Freund [33j, Kostrov g9l [50], and Willis [51J. It strictly applies to a 
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perfectly straight semi-infinite crack in an infinite plate propagating at a steady or non- 
steady rate with loads applied to the crack faces. The latter can can be mapped 
by superposition to a broader class of loading configurations, taking advantage of the 
linearity of the field equations and boundary conditions in ^ and ^ . This calculation 
is also valid for bodies of finite size for times that are shorter than the time needed 
for waves to bounce off of the boundaries and return to interact with the crack. With 
these restrictions the calculation is exact, and holds with remarkable generality. It is 
commonly thought that this is the only equation of motion for a dynamic crack, but 
this is not the case (see below). 

The full details of the derivation can be found in Freund's book [33] and with a less 
detailed review in [52J. The asymptotic displacement fields u for a crack propagating 
in an infinite plate have previously been presented in equation as a function of 
the stress intensity factor Kj(v). The dynamic theory predicts that the dynamic stress 
intensity factor Kj(v) can be written as 

K^v) = fc(v)if 5 (loading, I) , (13) 

with 

k(v) ~ (1 - v/c R )/y/l-v/c d . (14) 

k(v) is a universal function of the crack growth rate v — l alone. K s (loading, I) depends 
on the loading configuration and the crack length /, but not on v [33J. The computation 
of the latter might be difficult and may require numerical techniques, but it involves 



no dynamics. Equation (13) is remarkable because it decomposes a dynamic quantity 
Kj(v) into a purely dynamic universal function k(v) and a problem-specific quantity 
K s (loading, I) that does not involve dynamics at all (e.g under applied tensile stress 
far from the crack faces cj OQ) K s ~ a^Vl)- Thus, a dynamic problem is reduced to a 
non-dynamic one. 

For a crack moving at any velocity v in an infinite plate, the energy balance of 



equation (12) generalizes to 



1 - v 2 ft - u 2 )K 2 (l) 
I» = G(v) = -^-K*{v)A{v) = { - J->Uk?( v )A(v) , (15) 

where A{v) is yet another known universal function [33] and for simplicity we suppressed 
the loading dependence of K s . To a good approximation we have k 2 (v)A{y) ~ 1 — v/cr, 
which yields [33] 



2 



T(v)~^^K' 2 s (l){l-v/c R ) v~c R 



1 - 



Y{v)E 



:i-^)m) 



(16) 



The crack growth equation in (16), which is an equation for the time evolution of the 
crack tip location /(i), has the following rather interesting features: 

• The quantity G s (l) = (1 — v 2 )K 2 (l)/E can be interpreted as an effective 
thermodynamic force for crack motion. No motion takes place if G s < T. For 
r —> (no resistance to crack propagation) or G s — >■ oc (infinite effective driving 
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force), a crack will accelerate to a finite limiting speed equal to the Rayleigh wave 
speed cr. 



Equation (16) depends on the instantaneous tip location I and its instantaneous 
speed v — /, but not on higher order time derivatives. In particular, the absence 
of an inertia-like acceleration term v = I in this equation implies that the crack tip 
can be treated as a massless particle/defect that responds instantaneously to any 



change in either T or the driving force G s . From this perspective, equation (16) 
takes the form "v = F/rf\ where the RHS is the ratio of the net driving force a F" 
over an effective mobility u rf\ A term u mv" , where "m" is an effective mass, is 
missing from this equation. 



Comparison of equations (11) and (16) reveals that the dynamic contributions to 



the energy release rate serve to reduce its value relative to its quasi-static limit, 
G(v = 0). This is basically due to the increase in kinetic energy required to move 
the surrounding material away from the crack faces as the crack tip progresses. 
This energy cost diverges as v — >• cr and is the heart of the existence of a finite 
limiting velocity for a crack. 



1.5.2. The motion of a crack in an infinitely long strip When a crack is propagating 
within an infinitely long strip, its motion is described by a qualitatively different 
equation derived by Marder [47] . This crack configuration is commonly used to obtain 
model-independent measurements of T{y) [53j [54J [55]. The rationale behind this is as 
follows. Consider a semi-infinite crack moving at a constant steady-state velocity v in 
an infinitely long strip of width 2b (in the y direction) along its symmetry axis (in the 
x direction). Energy is stored in the strip by displacing its boundaries at y — ±6 by 
a constant amount. Far ahead of the crack tip (x —> +oc) the strip stores a constant 
energy per unit length, W ) whereas far behind the crack tip [x — >• — oo) all of this 
energy has been transformed into creating new crack surfaces. Under these conditions 
energy balance requires that Y{v) — W since translational invariance simply implies that 
when a unit length of new surface is formed at the crack tip an amount W of energy is 
released at x +oc while a new unit length of relieved stresses is added at x — >• — oo. 
Marder considered a crack's dynamics within the strip under non-steady conditions 
[47] . Performing a perturbative analysis where the dimensionless acceleration bvj(? d was 
assumed to be small, energy balance yielded 



G ~W 



1 



~ W 



1 



bv 



1 



= i» 



(17) 



Let us compare equations (16) and (17). One obvious difference is that the former 



depends on the crack's tip position /, while the latter is independent of it, but rather 



depends on the strip's half width b. More importantly, in contrast to equation (16), 



equation (17) tells us that the motion of a crack in a strip explicitly depends on its 
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acceleration v. Rewriting equation (17) as 



we observe that this equation takes the form "m(v)v = F", where the net force "F" 
is proportional to the difference between the effective thermodynamic force W and the 
crack propagation resistance T(y), and the velocity-dependent effective mass "m(v)" 
is proportional to f(v) ~ (1 — v 2 /c 2 R )~ 2 . One family of steady state solutions is 
obtained when T{v) = W . However, as T{v) is a material function that is bounded 
for < v < cr, one can externally set W such that W > T(v) for every < v < cr. In 
this case, the solution approaches v^cr in such a way that the effective mass diverges, 
m(v) ~f(v) — >-oo, and v^O while their product remains finite. 

While both equations (16) and ( fl7| ) (or equivalently equation (JT8]) ) are 



manifestations of energy balance at the crack tip and both predict v^Cr for sufficiently 
intense loadings, the underlying physics is quite different. In an infinite system elastic 
waves always outrun the crack tip whereas in a strip these waves are reflected back into 
the system and forced to interact with the crack. This qualitative difference in crack 



dynamics predicted by equation (16) and the dynamics that occur when a crack is able 
to interact with its "history" and acquire effective "inertia" (e.g. as in equation Jl8|)) 
is often overlooked. 

1.6. Using soft materials to test LEFM 



Whereas the LEFM predictions of fracture onset (e.g. equation (12)) have been 
experimentally validated for many materials (in fact this equation is used to measure 
the fracture energy Y at the onset of fracture), only a few direct experimental tests 
of the LEFM predictions for the dynamic behavior of cracks have been performed 
[561 157] . Historically, early measurements showing the inability of the maximal observed 
propagation velocity of a crack in brittle materials to approach cr had suggested that 
something was seriously wrong with our fundamental understanding of rapid crack 
dynamics, but it was unclear where the root of the problem lay. Later on, the problem 
was attributed to the existence of instabilities, which were observed both experimentally 

and numerically/theoretically [551 E2 ESI ESI 

The reader is referred to [52] for a comprehensive review of some of these works. 

The main experimental obstacle to direct tests of the theoretical predictions is the 
very high propagation velocity of cracks in typical brittle materials. In brittle polymers 
c s surpasses 1000 m/sec and in soda-lime glass (window-pane glass) shear velocities 
surpass 3000 m/sec. At these extreme velocities, direct visualization of displacements 
or material deformation in the vicinity of the crack tip is extremely difficult, necessitating 
extremely rapid high-speed cameras with frame rates in excess of 10 6 frames/sec. While 
such high speed cameras exist, measurement of near-tip crack deformations is further 
complicated by the microscopic size of the singular region in such materials. For 



The Dynamics of Rapid Fracture: Instabilities, Nonlinearities and Length Scales 13 



Table 1. 


Material properties of common 


brittle materials. 


Material 


Young's modulus (MPa) 


Rayleigh wave speed, cr (m/sec) 


Soda-Lime glass 


70000 


3400 


PMMA 


4000 


950-1200 


Polyacrylamide 


0.1 


5.6 



example, a sheet of brittle glass typically fractures at strains of less than 0.1% which in 
a sheet of 200 mm extension would result in a total crack tip opening of ^100— 200/xm. 
Moreover, the singular region near the crack tip will be of the order of 1-10 nm. Thus, 
real-time visualization of the singular part of the strain field surrounding the tip of a 
crack at reasonable spatial resolution would necessitate much better than 1/im spatial 
resolution at frame rates exceeding 10 6 . This is prohibitively difficult with existing 
imaging technology. 

The above obstacles can be circumvented by the use of soft materials (E « 100 
kPa), such as aqueous gels. These neo-Hookean materials are both elastic and can be 
made to be brittle, by tuning their chemical composition [63]. The shear wave speed 
of soft Polyacrylamide gels (13.8% acrylamide, 2.6% bis-acrylamide) can be tuned to 
be approximately 6 m/sec, nearly 3 orders of magnitude slower than "standard" brittle 
materials (see Table [I]). In addition, experiments using polyacrylamide gels indicate 
that these materials can be made to be relatively "tougher" than other brittle materials; 
typically fracturing at 5-10% strains [63], I94J [95] . Under these conditions, the near-tip 
singular region can exceed 1 mm in size thus enabling direct visualization by means of 
moderately fast video cameras. 

Can these "novel" materials indeed be used to investigate the brittle fracture 
of "standard" materials? The answer is, emphatically, yes! As figure [T] clearly 
demonstrates, crack dynamics are in excellent agreement with both the predictions 
for an infinite medium (equation (16)) and for propagation in an infinitely long strip 
(equation (JTtJ) ) . In the experiments on gels described in figure [TJ crack instabilities 
were suppressed (see section 2.2) so as to enable single cracks to attain velocities in 
excess of 0.9cr. Experimental details can be found in [95]. It is worth noting that 
the experiments in these comparisons were performed with no adjustable parameters. 
Excellent quantitative agreement was obtained for the crack dynamics in both cases as 
well as for the universal functions k(v) and f(v) that were predicted in, respectively, 
equations (16) and (17). Thus, not only do these types of materials correspond precisely 
to LEFM predictions for the motion of a single crack, these measurements have provided 
the most precise testing grounds over the widest range of crack velocities to date. 

The excellent compliance with LEFM predictions with experiments serves to 
validate both equations of crack growth rate (16) and (17), and could essentially serve 
as a "definition" of what we mean by a brittle material. The underlying assumptions of 
both equations of crack growth rate are: 
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Figure 1. Both equations of motion predicted by LEFM for an infinite medium 
(equation ( fl6| )) and for an infinite strip (equation (JTrj) ) are in excellent quantitative 
agreement with experiments performed using a polyacrylamide gel [55] . (a) 



Comparison of equation (16) for crack dynamics in an infinite medium for different 
applied values of a^. The inset shows a comparison of the universal function k(v) 



in equation (14) with experiments, (b) Crack dynamics for a finite strip compared to 
infinite medium prediction (solid line) as a function of the crack length, Z, normalized 
by the strip width, 2b. Once I w b the dynamics change dramatically and the crack 
dynamics are no longer consistent with equation (16) (red line), (inset) A comparison 
of the function f(v) predicted in equation (17) [47 (solid line) with experimental 
measurements in a strip geometry, (c) a sequence of photographs of the crack tip 
profile during the transition from an effectively infinite medium to a strip geometry. 
Note how the parabolic form of the tip at I ~ 0.66 (top) transitions at I ~ 1.46 (center) 
to a "tadpole-like" form at I ~ 2.2b (bottom) as the crack interacts with the waves 
reflected from the vertical boundaries. The parallel crack faces well behind the tip 
correspond to the displacement of the vertical boundaries. The arrow length is that of 
the half-width, 6, of the strip. See [55] for details of the experiments. 



• Only a single "simple" straight crack is propagating in the material. 

• Small-scale yielding is obeyed in the material, i.e. all of the dissipation that takes 
place within the fracture process zone is confined within the small singular region 
surrounding the crack tip. 

• The energy dissipated within the region per unit crack extension (the fracture 
energy T) is a material-dependent function of the instantaneous crack velocity v. 
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2. The limitations of LEFM: Crack instabilities 

2.1. The micro-branching instability 

In the previous section we demonstrated that, as long as fracture occurs via a single 
"simple" crack, we attain an excellent picture of its overall dynamics by balancing 
the energy driving it, G, with the dissipation at its tip, Y. Unfortunately, this is not 
a complete picture of brittle fracture. Experiments in a variety of different brittle 
materials have shown that beyond a critical velocity of about v c w OAcr a single 
crack becomes unstable via an instability coined the "micro-branching" instability 
[521 1581 1591 I6T1 l60l 1631 [96]. Beyond v c a single crack can undergo spatially local 
"frustrated" crack branching events, where a single "mother" crack gives birth to 
successive short-lived "daughter" cracks. The small "micro-branches" formed by the 
daughter cracks are confined to a very short range surrounding the crack tip (1— 100/im in 
PMMA). The micro-cracks become progressively longer and more dense with increasing 
G. In fact, sufficiently beyond v c the amount of net fracture surface formed by both the 
mother and daughter cracks can increase extensively; a fracture surface increase of over 
an order of magnitude was measured in PMMA as v increased from v c to 1.5v c [5H 155] . 



A brief review of this instability is presented in figure [2| where typical examples 
of fracture in both soft gels and glass are described. The figure clearly demonstrates 
the degree of universality of this instability. Figure [2^i shows that nearly every feature 
that characterizes the fracture process in glass is quantitatively similar (when properly 
normalized by the relevant velocity and length scales) to those features observed in 
the fracture of brittle gels [63] . These features include both the functional form of the 
subsurface profile of the micro-branches and the corresponding structure observed on the 
fracture surface that is formed once the instability is excited. In both materials, micro- 
branches are highly localized in the z direction (the depth of the sample, perpendicular 
to both the propagation direction x and the loading direction y). In addition, micro- 
branches form directed chains along the fracture surface in the direction of propagation; 
each micro-branch triggering the next one along the chain. This localized structure along 
z (crack front direction) highlights the intrinsically 3D nature of the micro-branching 
instability. This may explain why descriptions of the instability using generalized 
energy balance arguments for branched cracks in 2D media have not yielded quantitative 
agreement with measurements of instability thresholds [89], [97] . 

The slowing down of the fracture process that was so useful in quantitatively testing 
the predictions of LEFM for "simple" cracks (e.g. figure [j]) is also a distinct advantage 
in investigating the properties of this instability. A typical series of photographs of the 
instantaneous crack-tip profiles that bracket the onset of the micro-branching instability 
is presented in figure[2b. Prior to the instability onset, the crack tip propagates smoothly 



and is described well by the parabolic form predicted by equation 10, When micro- 



branching takes place, the overall propagation becomes jerky and the form of the tip 
becomes blunted, as the branches compete with the main crack, before being "outrun" 



The Dynamics of Rapid Fracture: Instabilities, Nonlinearities and Length Scales 16 





Figure 2. Universality of the micro-branching instability (a) A comparison between 
the instability in both soft polyacrylamide gels (left) and soda-lime glass (right). Once 
a crack attains a critical velocity of v ~ 0.4cr, a single "simple" crack may become 
unstable to the micro-branching instability. At this point (top) the instantaneous 
velocity of the crack undergoes violent oscillations that correspond to the formation of 
subsurface frustrated micro-branches (center), where side views (xy plane) of micro- 
cracks in both materials are presented (x is the propagation direction, y is the loading 
direction and z is the direction along the crack front). Micro-branches are very 
similar in both materials, characterized by a power-law functional form, (bottom) 
Photographs of the resulting fracture surface (xz plane) formed by cracks propagating 
from left to right at v ~ 0.5cr. The chains of structures on the fracture surface 
correspond to chains of micro-branches that are aligned in the propagation direction 
x and highly localized in the z direction. The vertical dimensions of the photographs 
are about 0.5mm. The cartoon in the center describes the geometry; the xz fracture 
surface is formed by a propagating cracks whereas micro-branches (dotted lines) extend 
below the surface, (b) A series of photographs of the quasi-parabolic profile of running 
cracks in polyacrylamide gels. The two lower left panels depict the parabolic crack 
tip opening profile of "simple" cracks propagating just prior to the onset of the micro- 
branching instability, as predicted by equation (10). The photographs correspond to 
the green dot on the velocity vs. crack length figure in the upper right corner. Right 
panels: 4 profiles of cracks undergoing the micro-branching instability, corresponding 
to the red dots in the velocity measurements. Note the small micro-cracks formed at 
the tip of the main crack. 
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by the winner. 

2.2. Oscillatory instability in rapid fracture 

The use of gels has enabled us to both observe and manipulate cracks at a level that 
is unprecedented in the study of fracture in "standard" brittle materials. For example, 
in gels it is possible to drive the system at very large strains, prior to fracture, and 
therefore achieve high acceleration rates. Applying these high accelerations, Livne et al. 
[63] found that, while the value of v c /cr at the lowest acceleration rates corresponded to 
the value observed in PMMA and glass, v c is - on average - a roughly linearly increasing 
function of the acceleration. This same study also established that the transition to 
micro-branching is highly hysteretic with features that are characteristics of an activated 
process accompanying a first order phase transition. A typical activation of a micro- 
branch chain is demonstrated in figure [3^. Beyond a minimum value of v c ^QAcr ) there 
is a bistable region of velocities in which either a single or multi-crack state can exist. 
In this region, the instability may be triggered when random perturbations surpass a 
critical threshold for activating the first micro-branch. This yields a finite probability 
to bifurcate in each time interval for v>0Acr. 



Figure 3. The micro-branching instability can be suppressed by using very thin 
samples and high accelerations, (a) Micro-branching is an activated process. The 
arrow indicates a chain of micro-branches generated at a nucleation site, (b) Once 
a chain of micro-branches encounters the free surface of the sample (arrow), it is 
annihilated. In both (a) and (b) fracture surfaces of width 0.5 mm, created by cracks 
propagating from left to right, are shown. 



The sub-critical nature of the instability provides us with a means to achieve single 
crack states at unprecedentedly high velocities [98] . The micro-branching instability 
is suppressed by driving cracks at very high acceleration rates. Furthermore, once the 
instability is activated, a chain of micro-branches that encounters a plate boundary (free 
surface of a plate) will disappear as shown in figure ^p. Thus, when an experiment is 
performed using very thin plates, activation centers that can trigger the instability are 



Z 




X 
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relatively sparse and, if the instability is triggered, the resulting chain of micro-branches 
quickly encounters a free surface at the edge of the fracture plane and disappears. This 
"trick" was used to obtain the unprecedented range of crack velocities used to validate 
the LEFM equations of motion in figure [TJ 

Despite these "tricks" , a close look at figure [T] indicates that it is still impossible to 
reach the asymptotic velocity of cr. When the micro-branching instability is suppressed, 
a new and unexpected oscillatory instability is observed at a critical velocity of about 
0.9c s [98|, [99]. Similar wavy cracks traveling at velocities beyond the shear- wave speeds 
have also been observed in experiments on latex sheets under biaxial tension at extremely 
high (> 100%) strains jlOOl 11011 1102j . The oscillatory instability in the gels, however, 
may be qualitatively different in nature. These (the gel) oscillations occur at clearly 
subsonic velocities and are driven by purely uniaxial tension at relatively small (~ 15%) 
strains. 



a 




Figure 4. When micro-branching is suppressed, a new oscillatory instability occurs 
at v = v osc = 0.9cr. (a) Successive xy profiles taken at 690/is intervals of a crack at the 
onset of the instability. The tip of an accelerating "simple" crack starts to oscillate 
in the y direction, forming a (b) (top) wavy sinusoidal pattern in the xy plane with a 
smooth and (bottom) featureless fracture surface. 

While an oscillatory instability has been predicted for rapid cracks J93J 1103] . the 
predicted wavelength in these models was expected to scale with the size of the system. 
Experimentally, none of the characteristic scales (e.g wavelength or amplitude) of the 
observed oscillations were seen to be at all dependent on sample geometries or dimensions 
[98] . It was suggested that this characteristic time/length indicates that a new intrinsic 
length scale was needed to describe these dynamics. This time/length scale, which 
could not be explained in the framework of LEFM, hinted that these new effects could 
be related to the process zone. 
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3. Beyond LEFM: regularization of the singularity and dissipation 

As we explained in detail above, LEFM is not a complete and self-contained theory of 
fracture dynamics. In particular it does not account for near crack tip/front dissipation, 
consequently it requires the fracture energy as an external input, and it does not account 
for crack path and velocity selection. There have been numerous attempts to address this 
problem by formulating models that go beyond LEFM in various aspects. The majority 
of this work has focused on purely 2D media, though not exclusively. A comprehensive 
review of these efforts goes well beyond the scope of this paper. For completeness, we 
briefly mention some of these works below. 

A central issue in the context of models of dynamic brittle fracture that go beyond 
LEFM is how to regularize the LEFM \j\fr singularity. Usually this regularization 
involves dissipation and hence it also gives rise to a fracture energy. A classical way to 
regularize the continuum LEFM singularity is to introduce a finite-size cohesive zone 
at the crack tip such that the singularity of the stress field is canceled out on the 
scale of the cohesive zone [341 1104j . Finite Elements Method (FEM) calculations have 
extended the cohesive zone idea to include cohesive forces between any two elements 
in the bulk. Cohesive elements calculations are very popular in engineering contexts, 
and have successfully reproduced several aspects of the dynamics of brittle fracture 
[SZl EHl EES EES UnZl EH EH HHJJ. Despite these successes, many aspects of these 
models are largely phenomenological. These include the introduction of external length 
scales near the crack tip. In addition, some inherent problems with such models are 
known to exist - see j3U H3 llllj . 

Much insight has been obtained by the use of lattice models with simple interaction 
laws, where the lattice constant provides a regularization length scale for the LEFM 
singularity [5llMIMl7^l85i rm[T^ 

11241 11251 1126j . It is difficult, however, to directly relate these models to the failure of real 
non- crystalline solids, where a natural regularization scale is not obviously apparent and 
well-defined slip systems do not exist. Recently, a continuous random network model 
has been developed and studied in detail j!27j . Direct molecular dynamics simulations 
have also successfully reproduced some of the phenomenology of dynamic brittle fracture 
ESI EE 11281 11291 113()j . Like lattice models, these calculations are generally performed 
in an ordered atomic material using largely phenomenological potentials. Indeed, 
attempts to use realistic crystalline potentials (e.g. in the fracture of silicon) have 
failed to quantitatively reproduce experiments jH[26]. Later on, it has been shown that 
quantum mechanical contributions to these potentials can not be ignored [271 EH] . 

Diffuse interface phase-field models provide a methodology to regularize the near- 
tip singularity. These models are self-consistent continuum formulations of brittle crack 
problems that incorporate an auxiliary phenomenological "phase field" to link the near- 
tip behavior to the linear elastic fields far from the tip |131l 11321 11331 1134] . These 
models provide a useful mathematical machinery that allows a self-consistent selection 
of the crack's speed, tip/front shape and path [ml ESI ESI ES EH E23 EE]. In 
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particular, they introduce a regularization length scale for the LEFM singularity and a 
time scale for near-tip dissipation. Phase transformation models of fracture were also 
studied in detail |141i I142j . While largely phenomenological at present, this class of 
models is promising and has been useful in elucidating several important aspects of 
crack dynamics. These include laws of crack motion in isotropic and anisotropic quasi- 
static 2D fracture [1431 1 144] and the evolution of complex quasi-static crack patterns in 
3D under mixed-mode loading |145i I146j . 

Many of the models that go beyond LEFM have addressed the issue of crack 
branching. As we stressed above, the vast majority of these focused on 2D where 
a crack tip bifurcates into two tips ("macro-branching"). These include lattice models 
[721 18TI I112j , a continuous random network model j!27j , molecular dynamics simulations 
[66], phase-field models [12JJ H35l H36J HSH E38l E39J ED], phase transformation 
models [1411 1142] and cohesive zone formulations coupled with finite element calculations 
[69; 109 . II 10] . The intrinsic 3D nature of the instability, however, has until now escaped 
our theoretical understanding |145U146U l47j. While a number of explanations for the 2D 
instability have been proposed [701 EH [86], [89l [971 [141J , the precise physical mechanism 
leading to this instability has not yet been established. 

Finally, the role of elastic nonlinearities in dynamic fracture has also been the 
subject of recent research [TDJ [72 EE EB EH [US EHl USS ESQ E52] . The basic idea is 
that the LEFM singularity implies the concentration of large deformations near the crack 
tip that must invalidate the assumption of a linear elastic behavior. On quite general 
grounds one expects linear elasticity to break down in favor of nonlinear elasticity before 
irreversible processes set in. Based on this idea, a weakly nonlinear elastic theory of 
dynamic fracture has been developed. This will be discussed in detail in the next section. 

4. A weakly nonlinear theory of fracture: The effects of nonlinear elasticity 
near the crack tip 

In section [2| we described two different instabilities that arise in brittle fracture. Neither 
of these can be understood in the framework of the fracture dynamics that are embodied 
in LEFM. We have also shown that, as long as a "simple" crack remains stable, crack 
dynamics are in excellent quantitative agreement with LEFM, as demonstrated in figure 

What is going on? Energy balance is at the heart of the single crack dynamics. 
LEFM uses an assumption of linear elasticity (nearly) everywhere to calculate the flux of 
energy transported from the external loading by means of the linear elastic fields (where 
materials do not undergo large deformations) to the region of a crack's tip. On the one 
hand, the highly deformed regions that exist in the near vicinity of a crack's tip can not 
be dealt with in the framework of LEFM. On the other hand, as long as these highly 
deformed regions are "passive" and exhibit no important dynamics of their own, they 
can be incorporated into the assumption of small-scale yielding and crack dynamics can 
be completely described in the LEFM framework. The unexplained "intrinsic" scale that 
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lies at the heart of the oscillatory instability, described in the previous section, however, 
provided a hint that something might be fundamentally missing in this assumption. 
An equation of motion for a crack should determine both the crack growth rate and 
its direction of propagation [82], 11441 11531 1154] . It is crucial to understand that energy 
balance, which is behind the crack growth rate equations (16) and (17), provides us 
with no information whatsoever about the direction of crack propagation, only about 
the speed of propagation. Without an equation that determines a crack's path, the 
question of path stability cannot even be formulated mathematically. Intrinsic length 
scales might be at the heart of such an equation. 

The Kjyjr scaling of the stress field predicted by LEFM as the crack tip is 
approached is a power-law, which possesses no intrinsic scale. That is, our fundamental 
understanding of fracture is based on the assumption that fracture dynamics are 
governed by the singular stresses that dominate all contributions to the stress field 
up to the ("single") point where the material is broken. All of the details of how a given 
system is loaded and what the geometry of the sample and the crack is, are incorporated 
into a single scalar quantity, the stress intensity factor, K. This is sometimes called the 
assumption of "autonomy" , which means that the mechanical state within the near-tip 
region is uniquely determined by the value of K and is otherwise independent of the 
applied loadings and the geometric configuration in a given problem [43j. In particular, 
the length scale inherited in K (whose dimension is stress times square root of a length) 
is extrinsic, associated with either the geometry or the loading of the problem. Were 
this framework capable of explaining the oscillatory instability discussed above, it would 
have inevitably implied that the wavelength of oscillations is determined by an extrinsic 
length scale. As the experiments clearly show that this is not the case, we must therefore 
look for new physics and a new length scale in the region where the assumption of linear 
elasticity breaks down, i.e. in the near-tip region. 

The use of soft gels provides us with a unique experimental opportunity to explore 
this elusive region. The slow propagation velocities of dynamic cracks within these 
materials have enabled us to perform measurements of unprecedented accuracy of the 
material deformations in the near-tip vicinity of truly dynamic cracks. This allows 
us to quantitatively examine the form of the deformation fields at scales surrounding 
the crack tip where the assumption of linear elastic behavior breaks down. Before we 
discuss in detail the outcome of these experimental investigations, we will first describe 
a new theoretical framework that goes beyond LEFM. This will set the stage for the 
quantitative analysis to follow. 

As our goal is to go beyond LEFM, we should ask ourselves what is the first physical 
process that intervenes when LEFM breaks down near a crack's tip. LEFM is based on 
the assumption of linear reversible (elastic) deformation. While traditional approaches 
focus on irreversible deformation near crack tips, e.g. plastic deformation [15511156] . our 
premise is that first linearity breaks down, while the deformation remains reversible. In 
a nutshell, we maintain that the harmonic (quadratic) approximation for the deviation 
from a stress-free configuration - which is at the heart of LEFM as in equation ^ - 
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must first give way to anharmonic (nonlinear) corrections that do not necessarily involve 
irreversible processes (e.g. particles rearrangements, decohesion etc.). Below we explore 
the theoretical implications of this idea. 

4-1- Finite deformations: a tutorial 

To set the stage for theoretical developments based on nonlinear elasticity, we will 
need some background in the theory of finite elastic deformations. This theory is well- 
developed and is summarized in various textbooks, see for example [1571 I158j . To 
render this paper as self-contained as possible, we briefly review some basic elements of 
this theory. We start with equations @-([2]) and define the deformation gradient F as 
F = V<fi = I+H. The Green-Lagrange tensor E is defined as 

E = h^H + H T + H T H) = e + \h t H , (19) 

where the linear strain tensor e is given in equation ([3]). E is a rotationally invariant 
tensor that measures the change in length of material elements and is evidently nonlinear 
in H. The nonlinear elastic energy density functional £/*, i.e. the elastic energy per unit 
reference volume, of isotropic materials can be expressed as 

U = U(I e JIeJIIe) , (20) 

where {Ie^He^IHe} is a se t of invariants which we take here to be the principal 
moments 

I E = trE, II E = trE 2 , III E = trE 3 . (21) 

The first Piola-Kirchhoff stress tensor s, that is work-conjugate to the deformation 
gradient F, is defined as 

s = = F . (22) 

OF dE K } 

s quantifies forces in the deformed configuration per unit area in the reference 

configuration. As will be shown below, this quantity is very useful in calculations. The 

Cauchy stress tensor <r, which is the true mechanical stress that appears in momentum 

balance equations, can be expressed in terms of s as 

sF T 

a = . (23) 

detF v ; 

The linear momentum balance equation is formulated in the deformed configuration in 
terms of the Cauchy stress <r, cf. equation The differential operators are understood 
to be defined with respect to the deformed coordinates x r , even though in the context 
of linear elasticity in which these equations were introduced, we made no distinction 
between the reference and deformed coordinates. The momentum balance equation 
can be rewritten in terms of the first Piola-Kirchhoff stress tensor s and the reference 
(undeformed) coordinates x as 

\7 X • s(x, t) = po(x)d tt u(x, t) , (24) 
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where po(x) is the (time-independent) mass density in the reference configuration. 
Angular momentum balance, <t = <t t , can be expressed as 



sF T = Fs T . (25) 



It is important to note that equation (24) is defined with respect to a fixed, known 



domain. This is a great simplification when compared to equation ^ , which is defined 
with respect to an evolving, yet unknown domain. This feature makes the first Piola- 



Kirchhoff stress tensor s of equation (22) very useful, especially for the application 



of these equations to fracture mechanics, where a crack introduces time-dependent 
traction-free boundary conditions described mathematically as moving branch-cuts in 



the undeformed configuration. Equation (24) will serve as a starting point for deriving 
the weakly nonlinear theory. 

4-2. The weakly nonlinear theory of dynamic fracture 

The weakly nonlinear theory of dynamic fracture incorporates (weak) elastic nonlinearity 
by expanding the equations of motion describing the deformation in a medium containing 
a propagating crack up to second order in the displacement-gradient tensor if, which 
is regarded as the small parameter in the problem. The starting point in the derivation 
is to write down the elastic energy density U up to 0(E 3 ) as 

U(E) = -A (trE) 2 + /1 trE 2 + ft (trE) 3 + f3 2 trE trE 2 + fotrE 3 + G(E A ) , (26) 

where A and /1 are the standard Lame constants and {ft, ft, ft} are the second order 
elastic constants. The second order elastic constants are basic physical quantities 
that represent the leading anharmonic contributions to the interatomic interaction 
potential. These anharmonic contributions are known to be the origin of many 
important physical properties and effects such as the Griineisen parameters, deviations 
from the Dulong-Petit law at high temperatures, thermal expansion and the existence of 
thermal resistance; their implications to dynamic fracture are discussed in this review, 
{ft, ft, ft} are, therefore, not free parameters, but rather quantities that are either 
calculated from a fully nonlinear elastic energy functional, if known, or measured directly 



in experiments (see section [4~3| below) . Equation (26) reduces to the linear elastic energy 
density equation Q under two conditions; the strain measure in equation (19) should 
be linearized with respect to H and the second order elastic constants should be set to 
zero. 

Like the first order elastic constants, where the pair {A, /1} can be replaced by 
equivalent pairs (e.g. the Young's and bulk moduli), the choice of second order elastic 
constants is not unique. For example, the Murnaghan coefficients {£, m, n} constitute 
an alternative set of second order constants |158j . They are simply related to {ft, ft, ft} 
defined above through 

^ = 3ft + ft, m = ft + 3ft/2, n = 3ft . (27) 



The Dynamics of Rapid Fracture: Instabilities, Nonlinearities and Length Scales 24 



It is important to note that U{E) in equation (26) is valid for a general isotropic 



material in 3D, where E is regarded as 3D tensor. In many cases, however, one is 
interested in situations in which the deformation state is 2D due to some simplifying 
physical conditions/assumptions. The only change to the above formalism would be in 



redefining the second order elastic constants {/3i, P2, Ps} i n equation (26). For example, 
under 2D plane-strain conditions [45j we have 



E 



^ El XX E X y 0^ 

, V 



(28) 



for which equation (26) becomes 



U 



2D 



{trE 2D f + ix trE\ D + (trE 2D f + 



2/5 2 



trE\ D , (29) 



where E 2 d is a 2D tensor composed of the non-vanishing elements of E in equation 



(28). To keep things as general as possible, we use the general 3D form of equation 
(26) to derive general results (i.e. not for a specific material or geometry) below, even 



when dealing with 2D deformation conditions. This simply means that when applying 
these results to a specific 2D problem, the second order elastic constants { /?i , /5 2 , /?3 } 
should be properly selected. For example, in the case of the plane-strain conditions of 
equation (29), we define 1 = p 1 +/?2/3, /?2 = 0, /3 3 = /5 3 + 2/3 2 /3} and use these instead of 
{ /?i , /3 2 , /?3 } in all of the results to follow. Another example, for plane-stress conditions, 
will be discussed in section [U 



To proceed, we use equation (26) in equation (22) and expand the result up to 
second order in H to obtain 

dU 



8 = (I + H) 



dE 



Xtrel + 2jie 



(30) 



-A tr{H T H)I + fji H T H + A treH + 2fiHe + 
3/?i {trefl + f3 2 tre 2 I + 2/? 2 tre e + 3/5 3 e 2 + 0(H 3 ) = s lst + s 



Jlnd 



Here s lst stands for the part of s that is linear in H and s 2nd stands for the part 
of s that is quadratic in H. Related second order stress tensors appeared in the 
literature, see for example |158i 11591 1160j . It is straightforward to show that the last 



expression automatically satisfies the angular momentum balance in equation (25) to 
0(H 2 ). Then, s of equation (30) can be used in the linear momentum balance in 



equation (24) to obtain the most general equations of motion for the displacement field 
u, consistent to second order in H. 



Equation (30) reveals the structure and origin of second order nonlinearities. The 



first two terms are linear in H and correspond to the usual Hooke's law for small 
deformations. The remaining terms are second order in H. These terms can be classified 
into two different types. The first type corresponds to second order nonlinearities 
that emerge from the fact that the basic rotationally invariant strain measure E is 
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quadratic in cf. equation (19). These second order nonlinearities do not involve new 



constitutive parameters (hence sometimes termed "geometric nonlinearities") and can 



be readily identified as the nonlinear terms in equation (30) that are proportional to the 



Lame constants A and \±. The remaining nonlinear terms in equation (30) emerge from 



the fact that the constitutive relation is nonlinear (hence sometimes termed "constitutive 
nonlinearities") and therefore are proportional to the second order elastic constants 
{Pii $2) fa}- I n the most general situation, we expect both types of nonlinearities to be 
relevant. 

We are now ready to derive the most general weakly nonlinear equations for the 
displacement field u and to apply these to an asymptotic expansion near the edge of a 
crack, where displacement-gradients become too large for the linear approximation to 
be accurate. To this end, we introduce a controlled expansion of u of the form 



0(e 3 ) = u 



(i) 



u 



(2) 



(31) 



where e is a measure of the magnitude of displacement-gradients, not to be confused 
with the small-deformation strain tensor e defined in equation ([3]). The expansion in 



equation (31) can be now substituted into the expression for s in equation (30) and 



then the equations of motion (24) can be expanded in orders of e. To first order in e we 



obtain the standard Lame equation, cf. equation ([6]). 
The second order equation in e takes the form 



r(i)l 



/xV V j + (A + /i)V(V • u^ 2) ) + T\u 
where T is a functional whose components are given by 

^[ti (1) ] = ^[ti«] + vS d [« (1)1 



Podttu 



(2) 



(32) 



(33) 



Equation (32) has the form of the Lame equation (equation (ph) with an added effective 



body loading given by JFJW 1 )] . The expression for s 2rui 5 defined in equation (30), is too 
lengthy to be presented here. Note that we already focus on 2D her^|} 



To apply equation (33) to the dynamics of cracks, we recall that a crack is 



defined as two surfaces that cannot support stresses and therefore is characterized 
by tract ion- free boundary conditions on its faces a^hj = 0, where ft' is the outward 
normal on the actual crack faces, i.e. in the deformed configuration. These boundary 
conditions can be rewritten in the undeformed configuration in terms of the first Piola- 
Kirchhoff stress tensor s as s ' ^ jf rt j = 0, where ft is the normal to the crack faces in 
the reference configuration. The latter introduces a great technical simplification as 
they are defined with respect to a known domain defined by the crack faces in the 



\ In the spirit of equation (31), the stress s can be also expanded in powers of e according to 
s ~es^ +e 2 s( 2 ) J rO(e 3 ) = s^ + s^ 2 \ Therefore, equation (32) is in fact of the form dj = pod u u[ 2 \ 
It is important to note the difference between s 2nd defined in equation (30) and s^ 2 \ The former is the 



second order contribution to s in terms of H, while the latter is the second order contribution to s in 
terms of e. 
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undeformed configuration. Note that since we will not be interested below in solving 
a global boundary value problem, but rather focus on a near crack edge asymptotic 
analysis, we do not explicitly consider the external boundary conditions imposed on the 
sample; the latter are needed to render the global boundary value problem well posed. 

The boundary conditions Sijhj = can be expressed explicitly for a propagating 
crack as 

s xy (r, Lp ±7r) s yy (r, Lp ±7r) . (34) 

The boundary conditions of equation ( |34| ), considered to first order in e, yield equation 
([8]). To second order in e, these boundary conditions read 

- /ir- 1 ^ - ,id r uf - S x [u^\ = , 

- (A + 2/x)r- 1 d^4 2) " A d r u x ] ~ S y [u^] = , (35) 

for Lp = ±7r[§} S[u^] is an effective surface force vector, which is quadratic in and 
whose components are given by 

<s> (1) ] = -<W=±0, 

Sy[u^]= -^<W = ±7T), (36) 



where s 2nd is defined in equation (30). 

Our goal here is to consider the order e 2 problem for the mode I (tensile) symmetry 
of fracture. The theory for mode II (shear) symmetry is presented in detail in [152J. 
The second order problem is defined by equations (32) and (35). Using the leading 
terms (proportional to Kj) in the first order displacement field of equations M3| in 



equations (|33j) and (36), we obtain 
S x {r, ±7r) = , 

IS? • (37 » 

where g{jp\v) and k(v) have the dimension of stress and the numerical factor 327T is 
introduced in order to be consistent with previously used definitions |149i 11501 [161J. 
k(v) is given explicitly as 

_ 16a 2 d v 4 (X + fi) 
K[V} - ciD\v) ' 

which happens to be independent of the second order elastic constants {^1,^2^3}- 
The vector function g((f',v) is too long to be reported, but it possesses the following 
symmetry properties g x (Lp]v)=g x (—Lp]v) and g y { ( p]v) = —g y {—Lp]v). We emphasize that 
these functions depend on material properties through the elastic constants, though we 
do not explicitly write this dependence in the functions' arguments. 

§ Note that in the spirit of the previous footnote, these boundary conditions are equivalent to 

(2) (2) 
Sxy (r, ip = ±7r) = and Sy y (r, (p = ±7r) = 0, respectively. 
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Together with the vanishing effective surface force component S x and S y (r, n) = 
<S y (r, _ tt) in equation (37), we see that the second order problem corresponds to a mode 
I linear elastic crack problem with an effective body force ^(r, ip) that scales as r~ 2 and 
an effective surface force S y (r, ±ir) that scales as r _1 , both are quadratic in Kj. 



The solution of equation (32) was derived in detail in [1491 1150] and takes the form 



u 



(2) 



u. 



(2) 



32tt// 2 



Alogr + — log 



9 • 2 

it sin 4 



+ 5a s logr + ^^logll 



sinV 



[-Aa d <Pd - + v)] 



327T/I 2 

The displacement-gradients derived from this solution 

Vu (2) ocl/r , 



(39) 



(40) 



all exhibit a singularity that is stronger than the 1 /yr singularity of LEFM. In addition, 
this solution is characterized by logr displacement contribution. Both features were 
directly confirmed in experiments |149l 11501 1161 j . Y((p;v) is an r-independent solution 



of equation (32) which does not satisfy the boundary conditions of equations (35). The 



Fourier series representation of Y((p;v) reads 



N 



T x (ip; v) ~ ^2 °n( v ) cos(n^) and Y y ((p; v) 



n=l 



N 

— ^2 d n (v) sm(mp) 

n=l 



(41) 



Specific solutions demonstrated that a small TV is sufficient to obtain accurate 
representations even at very high crack propagation speeds [149J. 



As Y((p]v) does not satisfy the boundary conditions of equations (35), it does 
not constitute a complete solution of the second order problem. The remaining part, 



cf. equation (39), comes from a solution of the homogeneous counterpart of equation 



(32) (i.e. the standard Lame equation) that is characterized by a 1/r displacement- 



gradients singularity. This property is needed in order to be able to satisfy the second 



boundary condition in (35), which contains an effective surface force that scales as 1/r. 



Substituting equation (39) into this boundary condition results in the following relation 
between A and B 

_ 2fiBa s - (A + 2/x)d^T y (7r; v) - k{v) 

A-(A + 2 M )a2 ' W 

This relation reduces to the one reported in |149j . where A = 2/x and Cd = 2c s were used 
for a neo-Hookean material, and k(v) was measured in units of ji. 



The parameter B seems to remain undetermined in the solution in equation (39). 



If true, this result is remarkable as it suggests that the concept of the autonomy of 
the near crack tip nonlinear region (see the discussion at the beginning of section [4] 
and [33j [43]) is not always valid. However, autonomy is a central concept in fracture 
mechanics [33, [43] and it is difficult to see how it could be violated. How, then, 



can one reconcile the fact that equation (39), with equation (42), satisfies the second 



order asymptotic boundary- value problem for all £>, but still be consistent with the 
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concept of autonomy? The resolution to this apparent paradox was discussed in detail 
in [1501 1161 j . It was shown that the 1/r singularity is unique in the sense that it 
gives rise to a spurious resultant (integrated) force in the crack's parallel direction 
(where no boundary conditions are imposed) and hence the resultant force balance is 
not automatically satisfied in this direction, even though the solution itself satisfies 
the asymptotic boundary- value problem. Since such an unbalanced/spurious force is 
unphysical (as no physical process, except for inertial motion, can produce such a net 
force), it must be eliminated by demanding that the resultant (integrated) Newton's 
equation (both per unit sample thickness) 

/TT f*T P7T 

SijUjrdLp = v 2 p / fdf / d xx Uidip = p { (43) 
-TT JO J-7T 

is satisfied. Here / is the net force per unit sample thickness acting on a line of radius 
r encircling a crack's tip and p is the time rate of change of momentum per unit sample 
thickness of the material enclosed in the circle. Note that the steady state relation 



dt = —vd x was used. The y-component of equation (43) is automatically satisfied due 
to the mode I symmetry. The x-component imposes a real physical constraint. It was 
further shown in |161j that unlike the \j\fr fields of LEFM (derived from equations 



£)])), the 1/r singular fields carry no net momentum rate, 

pf = v 2 p [ fdf [ d xx ufdip = . (44) 

J0 J-7T 



This result implies that 



/i 2) = 8$ nj rdtp = 0. (45) 



The latter constraint allows the determination of B. Therefore, B is uniquely determined 
once the LEFM asymptotic fields are known and the autonomy of the near tip region is 
retained. An example of calculating B was given in |161j . 

Before we turn to discussing experimental tests of this theory, we discuss some 
important issues related to it. 

4-3. The determination of the second order elastic constants 

An essential input to the weakly nonlinear theory is the set of second order elastic 



constants /3 2 , /3s} that were defined in equation (26). As discussed above, the second 
order elastic constants are directly related to the leading anharmonic contributions to 
the interatomic interaction potential that are known to be the origin of many important 
physical effects. We include this section since we believe that the topic of second order 
elastic constants is not well-known to the general reader. To both rectify this and 
to emphasize the fact that these are real physical constants that lead to real physical 
effects, we will now briefly review a number of methods that have been successfully used 
to measure these quantities. 
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Methods for measuring the second order elastic constants are well-developed and 
rather widely used. Nevertheless, and in contrast to the first order elastic constants A 
and /i, their measured values for many materials are not well-documented and hence 
usually are not easily accessible. Here we briefly review some of the experimental and 
theoretical methods used for obtaining the second order constants. 

A direct theoretical method for determining the second order elastic constants 
becomes available once first-principles calculations of the structure and the interatomic 
interaction potentials of a given material are known. In that case, the elastic energy 
functional of equation (26) can be directly calculated and the elastic constants are 
obtained from a polynomial fit to the calculated energy-strain relation. In fact, yet 
higher order elastic constants can be obtained using this procedure. Recent first- 
principles density-functional theory calculations explicitly demonstrated this approach 
for single crystals [1621 1163] . Additional recent work employed tight-binding atomistic 
simulations to calculate the second order elastic constants of monolayer graphene 
j!64i I165J . Naturally, such first-principles approaches are more adequate for crystals 
with a well-defined symmetry than for amorphous solids. 

Some of the experimental methods for determining the second order elastic 
constants are described in the 1981 review paper of Hiki [166J, where a list of papers in 
which available data (as of 1981) is compiled. Static methods mainly involve measuring 
the onset of nonlinear variations of stress vs. deformation. Such static measurements 
for many materials have already been carried out in |167j . In addition to the static 
methods, dynamic methods exist that are mainly based on wave propagation. One 
such method [166J 11681 1169j is based on the idea that in a stress-free material, finite 
amplitude waves of a given fundamental frequency generate higher order harmonics 
(waves whose frequencies are integer multiples of the fundamental frequency) due to 
elastic nonlinearities. Second order elastic nonlinearities give rise to the generation of 
second harmonics whose amplitude depends on the second order elastic constants. This 
method has been applied to a number of materials, see for example [166] 11681 [169J. 

The most widely used method for measuring the second order elastic constants is 
acoustoelasticity, which is the acoustical analog of photoelasticity in optics. The method 
is based on measuring the speed of small amplitude plane waves that are superimposed 
on an applied static stress. The first expressions for the wave speeds as a function 
of the applied stress and the resulting second order elastic constants appeared in the 
pioneering work of Hughes and Kelly j!70j . Additional acoustoelastic techniques are 
discussed by Crecraft |171j . who developed a rather accurate acoustoelastic technique 
based on ultrasonic waves. 

Acoustoelasticity has been extensively applied to many materials, employing 
various techniques for measuring the wave speeds. For example, a Brillouin spectroscopy 
method was developed and applied to measure the second order elastic coefficients of 
solid polymers (172] . soda-lime-silica glass |173j and a bulk metallic glass |174j . Recently, 
a coda wave interferometry method was used to determine the second order elastic 
constants of complex solids such as concrete |175j . 
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Finally, the second order elastic constants can be determined by expanding a well- 
established fully nonlinear elastic energy functional. Such an energy functional usually 
becomes available by directly fitting experimental data for highly compliant materials, 
e.g. rubber-like materials and elastomer gels. In some cases the form of the energy 
functional is supported by a microscopic model, e.g. the neo-Hookean model, but in 
other cases symmetry considerations and pure phenomenology are sufficient, e.g. the 
Blatz-Ko energy functional that was used to fit some experimental data for foam rubber 
j!76i I177j . Once a fully nonlinear elastic energy functional is available, expansion to 
third order in the Green-Lagrange strain tensor E and a direct comparison to equation 
^ yields {Pufafo}. 

In Table [2] we present the second order elastic constants for several materials (a 
glassy polymer, two metals and a soft material) using available experimental data or an 
expansion of a known fully nonlinear strain energy functiona{jj] Table [2] reveals that the 



Table 2. The second order elastic constants /?2, ^3}, in units of the shear modulus 
/i, for various materials. The data for Polystyrene |170] , copper and aluminium |171j 
were obtained from acoustoelastic measurements. The constants for foam rubber were 
calculated using the fully nonlinear Blatz-Ko energy functional |176[ 1177] . 



Material Pi/v 

Polystyrene -2.56 -6.01 -2.41 

Copper 5.17 -3.73 -2.91 

Aluminium 2.29 -8.76 -3.32 

Foam rubber 1/6 -1 -8/3 



second order elastic constants have no definite sign and that they can be significantly 
larger than the linear elastic ones. While for (soft) foam rubber the second order elastic 
constants are of the order of the shear modulus, for "hard" solids they can be nearly an 
order of magnitude larger. Once /?2, ^3} are determined, the weakly nonlinear theory 
of dynamic fracture provides quantitative predictions that can be tested experimentally 
as described in section [H 

The second order elastic constants for an incompressible neo-Hookean material 



under plane-stress conditions will be discussed below (cf. equation (55)). 



4-4- Properties of the weakly nonlinear solution 

Some properties of the weakly nonlinear solution are highlighted below. These properties 
represent non-trivial extensions of the linear solutions and, as we will show in section [5| 
can be tested experimentally. 

• Crack tip profiles in the weakly nonlinear theory As was discussed above, an 
important and experimentally accessible property, is the crack tip opening profile 

|| The experimental papers reported the values of the Murnaghan coefficient's [158 , which were 



transformed into { fi\ , P2 ? P3 } using equation (27). 
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(the so-called "crack tip opening displacement" - CTOD). The weakly nonlinear 
theory predicts corrections to the parabolic CTOD predicted by LEFM, cf. equation 



(10). Using equation ( |39[ ), we obtain: 
K 2 



(A + a s B) log(r) + c n {v) cos (nir) 



X2 log (r) + xs 



which immediately implies 

(j> x {r, ±tt) = -xi (</> y (r, ±tt) T X4) 2 + X2 log [xi (^(r, ±tt) T X4) 2 ] + X3 • (47) 

The latter is valid in the weakly nonlinear region and describes a parabolic form 
corrected by a logarithmic stretch. Note that by setting X2 = X3 = X4 = 0, equation 



(10) is recovered. 



The sub-leading weakly nonlinear solution as a leading effect The second order 



displacement fields in equation (39) are the leading order corrections to the first 
order fields in equations ^ when a crack's tip is approached from the linear elastic 
region. There might exist situations in which these sub-leading weakly nonlinear 
contributions become the dominant ones, when the first order terms happen to 
vanish or become small for some physical reason. One such example was discussed 
in |152j . where it was shown that during mode II (shear) crack propagation, weakly 
nonlinear contributions are of tensile nature and hence may result in crack tip 
opening and tensile stresses ahead of the tip. This is a leading effect as these 
quantities are identically zero for the Mode II problem in LEFM. This effect may 
be relevant for problems like frictional sliding and the super-shear transition in 
mode II propagation |178i 1179] . 

Another example was briefly discussed in [148J 1149] . where it was noted that the 
tensile strain component d y Uy changes sign at a finite velocity, say ahead of a 
mode I crack's tip. Using equation ^ we obtain 

^■W = 0)= + (48) 

which can be easily shown to change sign from positive to negative at 

«b = \ Ljc\ + 8c2 - c\ . (49) 



For example, for q = 2c s we obtain v = 0.73c 5 , cf. figure |8p [148J. For v > 
LEFM predicts that d y Uy\r, (p = 0) is negative. This implies that, as the crack tip 
is approached, material points straddling y — come closer to one another instead of 
becoming increasingly separated as needed to precipitate fracture. This prediction 
of LEFM is not always appreciated. It was shown in [148] to be in contrast with 
experimental measurements in which d y u y y\r ) = 0) is always positive for v > vq. 
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As we will show in figure [8)3, the positive separation results from the dominant 
contribution provided by the weakly nonlinear theory; the second order contribution 
d y Uy\r,(p = Q) is positive and dominates d y Uy\r,(p = Q) in this velocity range. 
The length scale associated with the weakly nonlinear theory One of the most 
important aspects of the weakly nonlinear theory of dynamic fracture is that it 
introduces a new length scale into the fracture problem. This length scale will be 
shown below to be of prime importance for understanding (at least) one dynamic 
crack tip instability, and might hold the key for cracking other puzzles in fracture 
dynamics. The new length scale, which we denote by £ni(v), represents the scale in 
which the nonlinear elastic contribution to the mechanical fields become comparable 
to the linear elastic asymptotic contributions near the tip of a crack. Physically, it 
describes the scale at which LEFM breaks down. The subscript nl highlights the 
fact that the origin of this length scale is near tip nonlinearities. 
To understand the properties of the length scale £ n i(v) : we expand the deformation 
gradient as H ~ V x u^+V x u^ = H^+H^ near the tip of a crack. £ n i(v) emerges 
as a result of the different r-dependencies of and H ^ near the tip of a crack. 
According to equations ([9]), takes the following form 

H« = fc«( V /c a ,A/M)^= , (50) 

where h^(v/c S) A//i) is a calculable dimensionless tensorial function of the crack 
propagation speed v/c s and the ratio of the first order elastic constants A//i. 
According to equations (39), takes the following form 

JfW = h<?\v/c a ,\/ti, Pi/fi, fo/^fa/n)^- , (51) 

ji z r 

where h^ 2 \v/c S) A//x, ^//i, Ps/fj) is a complicated, yet calculable, dimension- 
less tensorial function of the crack propagation speed v/c s , the ratio of the first 
order elastic constants A//i and the dimensionless second order elastic constants 
{/3i/ /a, fa/ /a, fa/ /a}. £ n i{v) can be estimated according to 

\H(X =eni ~\HU\ r=eni , (52) 

which implies that scoling-wise it is given by 

k 2 r 

£ru~ ~i ~ - , (53) 

where we used the fact that the fracture energy T is proportional to K] / \i [33J [43] , 
(cf. equation (15)). It is important to note that l n \ is a dynamic length scale that 
evolves with the crack propagation speed v due to both the fracture energy T(v) 
and the ^-dependence of the pre-factor of equation (53). The pre-factor, itself, is a 
nontrivial function of both v and the first and second order elastic coefficients and 
is not necessarily of order unity. Incorporating l n \ into the theory of fracture, as 
will be discussed below, offers novel insight into the failure dynamics of solids and 
is one of the main take-home messages of this review. 
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5. Comparing the theory to direct experimental measurements 

We now wish to use the soft gels discussed above to directly test the predictions 
of the weakly nonlinear theory of dynamic fracture. To accomplish this, we must 
first determine the second order elastic constants { /?i , /3 2 , /?3 } for polyacrylamide gels. 
The neo-Hookean constitutive law provides a good description of the nonlinear elastic 
behavior of many compliant materials |180j . This constitutive law is a natural extension 
of Hooke's law to finite deformations [181J and is the simplest description of rubber- 
like behavior based on the Gaussian chain statistical model of entropic elasticity [180J. 
The gels used in the experiments described below are incompressible and deform under 
plane-stress conditions (thin samples). Under these conditions the neo-Hookean energy 
functional takes the form [182] 

U(F) = | [tr(F T F) + (det F)~ 2 - 3] , (54) 

where F is the 2D deformation gradient tensor. detF appears in the elastic energy 



functional due to the incompressibility condition. Equation (54) conforms with the 



general expression in equation (26) if one identifies A = 2/x and 



fa = -4/i/3, fa = 0, fa = -8/i/3 . (55) 

With these values of the fa's in hand, we can now compare the general weakly nonlinear 
theory directly to experimental measurements. Let us first consider the measurements 
of the CTOD in the near vicinity of the crack tip, as presented in figure [5} In figure 
[2)3 we saw that the crack tip can be described by a parabolic form, as predicted by 
equation (10). As shown in equation Q, the crack tip curvature is wholly determined 



by Kj and can be used to directly measure both Kj and the fracture energy T(v) (via 



equation (15)). Such measurements of T(v) are not simply an exercise in curve fitting, 
but, as shown in figure [5)3, are in excellent agreement with independent measurements 
of T(v) performed using a strip configuration [95] (cf. equation (jTTj) ) . 

Let us now look a bit more closely at the close vicinity of the crack tip. We find, 
as shown in the right panel of figure [5^i, that the CTOD predicted by LEFM fails to 
describe the crack tip within a region 5, defined as the distance between the real crack 
tip and the LEFM prediction. As LEFM predictions perfectly describe both the CTOD 
at intermediate scales (e.g. center panel of figure [5^) and T(v) (figure^), the existence 
of 5 is not an "artifact" of a "poor" fit to the data. In fact, figure [5J3 demonstrates that 
5(v) is a systematically increasing function of u, whose value increases by over an order 
of magnitude as v doubles. Comparison of the steep increase in S(v) with the mildly 
varying T{v) over the same velocity range, suggests that the increase in 5 is not due 
to an increased dissipation with i>, but is, instead, an effect of elastic nonlinearity. A 



large part of this effect, as described in weakly nonlinear solution, equation (47), is a 
"stretching" of the x axis due to the new logarithmic term resulting from weak nonlinear 
elasticity. This is demonstrated in figure |6| where we present a detailed comparison of 
LEFM predictions (dotted line) and those of the weakly nonlinear theory (red line) given 
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Figure 5. In the close vicinity of the crack tip the pure parabolic form of the crack 
tip opening displacement (CTOD) predicted by LEFM breaks down at a scale, S. (a) 
The crack tip opening at different scales: system scale (left) millimeter scale (center) 
sub- millimeter scale (right). The parabolic form predicted by equation 10 (dashed 
line) fails to describe the near-tip region at scales below 5, the distance from the real 
to the predicted tip of the crack, (b) The fracture energy, T, obtained both from the 
crack tip curvature (equation ([9| and using equation ( fl5| ) (circles) and from steady 
state measurements in a strip (squares) (see [95]). (c) S(v) for the same velocity range. 
These results suggest that S(v) is not a dissipative scale, as T(v) is fairly constant with 
v while S(v) increases by nearly over an order of magnitude over the same velocity 
range. 



by equation (47). For completeness, we show the concatenation of the weakly nonlinear 
theory to predictions of an asymptotic theory [94j I183j accounting for strong elastic 
nonlinearity (at scales within 100/xm from the tip) in neo-Hookean materials that are 
encountered at very large (3> 1) strains. The theoretical comparison described in figure 
[6^i was obtained from first principles using no adjustable parameters, as the second order 
elastic constants for neo-Hookean materials are known (cf. equation (55)). 

The same analysis that yielded figure [6^i for a particular crack velocity was 
performed for the range of v studied in [94]. This analysis resulted in the excellent 
quantitative agreement with the measured values of S(v) presented in figure^. Figure[6] 
therefore demonstrates that simply accounting for the nonlinear elasticity of the material 
surrounding a crack's tip provides us with a complete quantitative description of the 
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Figure 6. The nonlinear elastic contributions to the CTOD. (a) A comparison of 
measurements (circles) at v ~ 0.7cr of the near-tip CTOD with predictions of LEFM 



(dotted line), the weakly non-linear theory (equation (47)) (red line) and theoretical 
predictions for a neo-Hookean material (blue line) in the nonlinear "asymptotic" region, 
where strains are significantly greater than unity (see [94] for more details). Nonlinear 
elastic contributions to the CTOD provide a nearly perfect description of the form 
of a crack's tip, as long as crack instabilities are suppressed and a single crack is 
propagating, (b) The same nonlinear elastic theoretical analysis (solid line) provides 
an excellent quantitative description of measured values (squares) of S(v) [94] . 



CTOD of highly dynamic cracks. 

At what scale do elastic nonlinearities and dissipation take place in the neo-Hookean 
materials considered? Figure [7] describes direct measurements of the integrated energy 
flowing into contours of decreasing size that surround the tip of a crack propagating at 
v ~ 0.7c s [94] . The experimentally measured displacement fields enable us to directly 
calculate the energy flux through any contour C assuming a specific elastic energy 
functional [/, using the J-integral [HI [33] 

J = J i( U + lpodtUid t Ui) vn x + SijUjdtU^ dC , (56) 

where n here is an outward unit vector on C. The energy release rate is given by 
G = J/v [41] . This calculation puts a stringent constraint on U since only the physically 
correct one will make the result path (contour) independent, as must be the case as 
long as the path does not go through a dissipative region and the crack propagates 
at a steady velocity. The squares in figure [7] correspond to G for different contours 



obtained using the elastic energy function of equation (54), while the circles represent 
the corresponding calculations using the linear elastic approximation of equation Q. 
To the best of our knowledge, this is the first ever direct experimental estimate of the 
J-integral for a moving crack. This figure demonstrates three important things: (i) The 



elastic energy function of equation (54) properly describes the constitutive behaviour 



of this material down to the smallest scales near the tip of the crack since G is indeed 
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Figure 7. The energy flux G through different contours using the neo-Hookean 
energy functional in equation ( |54| ) (red squares) and its linear elastic approximation 
(blue circles). (Inset) The measured displacement field u x of a crack propagating at 
v~0.7c s (color bar in millimeters). Dashed blue rectangles mark every second contour 
used for calculating G. The value of G corresponds to the independently measured 
value of T(v c± 0.7c s ) (dashed black line). The results are presented in the reference 
(undeformed) configuration, where the crack is denoted by the black line with its tip 
at the origin. Adapted from [94] . 



shown to be path (contour) independent with a value exactly equal to the measured 
fracture energy T{v) at this propagation velocity, (ii) The progressive breakdown of 
LEFM as the crack tip is approached highlights the existence of the length scale, £ n i, 
associated with near tip elastic nonlinearities. At large scales (8-10mm) the deformation 
is small and nonlinearities in U make no difference (the circles overlap the squares). As 
the tip is approached (i.e. smaller and smaller contours are used), the deformation 
becomes progressively nonlinear and the LEFM predictions deviate significantly from 
the nonlinear ones (reaching about 50% deviation for the smallest contour of typical 
length of a few hundred fim. (iii) Down to the smallest contour, no dissipation is 
observed (G in the nonlinear calculation is still constant, cf. the squares). Hence, an 
upper bound for the dissipative scale in this material is a few hundred /xm. Comparison 
to the nonlinear "asymptotic" theory in further reduces the upper bound of the size of 
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the dissipative region to within ~ 20/xm from the crack tip 

We now wish to perform a direct comparison with the strain fields that are measured 
as close to the crack tip as possible. As we saw in equation (40), the nonlinear 
corrections to the strains have a stronger singularity (~ 1/r) than the ~1/ y/r singularity 
predicted by LEFM. In addition, at high velocities (~0.7cr) the nonlinear contribution 
to the strain becomes the dominant one (cf. equation ([48])) . In figure [8] we compare 
direct measurements of both u x (r,0) and e yy {r ) 0) = d y u y (r,0) to predictions of the 
weakly nonlinear theory (line) and LEFM (dotted line). The explicit form of the 
weakly nonlinear solution in equation (39) for a neo-Hookean material contains only 
two parameters (Kj, T) that cannot be determined from the asymptotic solution and 
therefore must be extracted from the experimental data. Kj was determined by fitting 
the far-field parabolic form of the CTOD to equation (46). The value of T for each v is 
unknown and was left as a free parameter [149J. We then follow [149J and use equation 
(39) with the nonlinear parameter B in equation (42) as a free fitting parameter to 
compare the theory to the measured experimental strains. For v = 0.20c s we also include 
the predicted curves using the theoretically calculated value of B. In figure [8] we present 
the resulting comparison for u x (r, 0) and e yy {r ) 0) with v/c s = 0.20, 0.53 and 0.78. 

It is evident from figure [8] that the agreement with the experimental data is 
excellent. For the low velocity curve there is essentially no difference between the results 
corresponding to the calculated and fitted values of B. The results clearly demonstrate 
the effect of the predicted \jr singular terms near the crack tip. In particular, the 
highest velocity (v = 0.78c 5 ) is larger than vo = 0.73c s calculated in equation (49), for 
which LEFM predicts that e yy {r, 0) becomes negative (cf. the dashed line in panel c). 
We see that the second order theory already avoids this intuitive conundrum for v>v . 
Thus, the second order nonlinear solution (solid line) both induces a qualitative change 
in the character of the strain and, moreover, yields excellent quantitative agreement to 
the measurements. This favorable comparison demonstrates that that the contributions 
of elastic nonlinearities are both important and generally unavoidable. This is especially 
true as high crack velocities are reached and the size of the nonlinear region, as evidenced 
by 5, becomes large. 



6. Cracking a dynamic crack instability 

In the previous section we demonstrated that the weakly nonlinear theory is a natural 
and a necessary extension of LEFM, providing an excellent quantitative description of 
the near-tip fields surrounding a propagating crack. These results provide a new and 
comprehensive picture of how remotely applied forces drive material failure in the most 
fundamental of fracture states: straight, rapidly moving cracks. In particular, these 
results reveal the beauty of how the hierarchy of linear and nonlinear elastic zones 
couple to transport energy from the macroscopic scales inherent in the external loading 
to the increasingly singular regions surrounding a crack's tip - before being dissipated at 
still smaller scales. On the other hand, as long as a single crack retains its path stability 
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Figure 8. Top: Measured u x (r, 0) (circles) fitted to the x component of equation ( [39] ) 
(solid line) for (a) v = 0.20c s with Kj = 1070Pav/m, T = -3150Pa and B = 18. (b) 
v 0.53c s with Kj = 1250Pav/m, T = -6200Pa and B = 7.3 and (c) v = 0.78c s with 
Ki = OSOPay^m, T = — 6900Pa and B = 26 Bottom: corresponding measurements of 
£ yy (r, 0) =d y Uy(r, 0) (circles) compared to the weakly nonlinear solution (solid lines) 
(cf. equation (31)) where T and B are taken from the fits of u x (r,0) and u y (r, n). 
(dashed lines) LEFM predictions (analysis as in |148j ) were added for comparison. In 
the left panels we added (dashed-dotted lines) the weakly nonlinear solution where 
the nonlinear parameters are theoretically calculated using Kj — 1040Pay / m and 
T = — 2800Pa. The results are nearly indistinguishable from those obtained by taking 
the nonlinear parameter B (in equation ([42])) as a free parameter (solid lines). Adapted 
from 



the theory provides nothing new with regard to a crack's dynamics; energy-balance 
essentially governs the overall behavior of a rapidly moving "simple" crack. 

Were dynamic cracks always simple cracks, we would now have a complete theory 
of fracture up to the dissipative scale. As we have shown in section [2j however, simple 
cracks become unstable in an number of ways. Furthermore, in the case of the oscillatory 
instability, a new length scale (the oscillations wavelength) was born, whose origin had 
nothing to do with all typical external scales in the system. We have shown that the 
weakly nonlinear theory provides a new and intrinsic length scale, t n \ (cf. equation (53)). 
This scale is determined dynamically; it represents the crossover length where nonlinear 
elastic effects become significant. Experimentally, t n \ is simply related to the scale 5(v), 
which describes the "extra" length that a crack tip receives due to the nonlinear elastic 
contributions. 

In a simple crack l n \ is a "passive" quantity. Although formed by the nonlinear 
elastic fields, it is simply carried along at the tip of the crack. In this section we explore 
the ramifications when this scale is no longer passive, becoming instead dynamic (or 
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active). When this occurs, non-trivial interactions can take place between the linear 
fields that drive the non-linear region and the non-linear response, as represented by 
t n \. We shall see that this active feedback between the linear and nonlinear regions 
provides the key towards understanding the oscillatory instability. 

6.1. A dynamic crack tip equation of motion and linear stability analysis 

Probably the first question one should raise in relation to dynamic fracture instabilities is 
why our understanding of them is much less developed in comparison to other, seemingly 
similar, instabilities in condensed matter physics and material science? For example, 
the micro-branching crack instability might appear somewhat similar to side-branching 
in dendritic crystal growth [185J. The linear regime of the solidification instability 
was essentially explained by Mullins and Sekerka nearly fifty years ago |185j . and that 
insight has been the basis for major advances in solidification theory and processing 
ever since [186J. This progress has stemmed directly from the fact that the dynamical 
evolution of the solid-liquid interface is governed on a continuum scale by a well-defined 
free-boundary problem. In contrast, it is not yet clear how to rigorously formulate an 
analogous free-boundary problem for fracture, where material failure is localized to a 
small singular region near the crack tip. Consequently we have, as of yet, no comparable 
understanding of dynamic fracture instabilities. One crucial missing ingredient is an 
understanding of the time and length scales associated with the physics near the tip of 
a propagating crack. In this section we will describe a recent attempt to derive a new 
dynamic crack tip equation of motion, incorporating the existence of the length scale 



The first step in this derivation closely follows the reasoning of |187j . Consider a 
crack in a 2D body, whose path is described by r tip {t) and whose tip is surrounded 
by a small nonlinear zone of scale i n \. We denote by t and n the tangent and normal 
unit vectors at the crack tip, respectively (see figure [9]). As we consider an arbitrary 
crack path under general external loading conditions, the LEFM near tip fields are 
generically characterized by both Kj (tension) and K n (shear) contributions. Consider 
then the discrete symmetry operation R n that transforms n —n. Under this symmetry 
operation, the relevant quantities of the asymptotic LEFM fields for both mode I and 
II fracture transform as follows: (i) Kj—^Ki (ii) Ku—t—Ku (hi) v—>v. 

Assuming material isotropy, one can write down the most general first order 
equations that are invariant under R n . The first equation is just a kinematic relation 
for the rate of crack growth 





(57) 



The second one describes crack tip rotation 



d t t oc K n (t)n 



(58) 



where the proportionality coefficient is a true scalar. 
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Figure 9. A crack with a small nonlinear zone of scale l n \. The angle 6 the crack 
makes with respect to the x-direction is shown at two times separated by a delay r^. 



The second step in the derivation amounts to estimating the proportionality 
coefficient in equation (58) by dimensional considerations. The existence of a length 
scale £ n i and a crack speed v suggests a time scale £ n i/v. Equation (12) allows us to 
define a typical stress intensity factor as K c 



(58) can be written as 



d t t 



' EY. Together, these imply that equation 
v K n {t) 



n 



(59) 



Li K c 

The adopted scaling approach assumes that all other material-specific properties of the 
nonlinear zone appear as a pre- factor of order unity in equation (59). Equation (59) 
can be rewritten in terms of the angle 6 that the unit tangent t makes with the x-axis 
as piUBB] 

v K n (t) 



d t 0(t) 



K r 



(60) 



The third step in the derivation follows from the observation that 0(t) is defined at 
the crack tip, while K n (t) is defined a distance ~ t n \ away from it. Therefore, causality 
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implies that Ku at time t cannot be affected by the crack faces created in the time 
interval [t— r^, t], with 

Td ~ Li/cni ■ (61) 
Here c n \ is the typical wave speed within the nonlinear zone, possibly of the order of the 
linear elastic wave speed c s , but not necessarily so. We note that using a single delay 
time Td is certainly a simplification of more complicated dynamics, but this simplified 
scaling assumption is expected to capture the essence of the physics involved. 

To formulate this idea precisely, we should express the physical K n at time t in 
terms of the mathematical JCi and Kn at a retarded time t—Td, taking into account the 
fact that the latter are defined with respect to a coordinate system rotated by #(£— t^), 
while the former with respect to a coordinates system rotated by 0(t), see figure [9} 
Kjit — Td) and JCn(t — T d ) are obtained from a pure LEFM problem with a crack path 
corresponding to t — Td and a simple consideration allows us to express K n (t) in terms 



of them. Equation (60), supplemented by the relation Knit) [/Cj(t— r^), Knit— Td)], 
constitutes the proposed dynamic equation of motion for the crack tip. It can be shown 
to reduce to the so-called "principle of local symmetry" - stating that cracks propagate 
so as to annihilate K u [US HEH1 EHS EHU EH E£3 E21 EEj - under quasi-static 
conditions. In this limit, the delay time Td is expected to play no important role. For 
propagation velocities of the order of the speed of information, however, new physical 



effects might emerge. It is important to stress that the derivation of equation (60) did 
not make any explicit reference to the origin of £ n i : which denotes the scale in which 
linear elasticity breaks down, either by nonlinear elasticity or by dissipative process. For 
example, in principle it could equally well apply to situations in which the breakdown 
of linear elasticity is dominated by plastic deformation. This remains to be verified 
experimentally in the future. 



Equation (60) is a nonlinear integro-differential equation for 9{t). With such an 
equation at hand, one can mathematically pose the question of the linear stability of 
the crack's path. To address this, consider a straight crack propagating at a steady 
velocity v under mode I symmetry conditions. Consider then a small perturbation of 
the straight path, characterized by an amplitude a and a wavelength A. In the limit 
a/A<l, we can consider linear modes of the form 6(t)~ae luJt /A, with $1(lj) = 2ttv/\, and 



linearize equation (60) with respect to a/ A. The resulting equation for uj determines the 
linear stability of the crack, employing the Willis-Movchan linear perturbation formalism 
[T3l [141 1154j . In particular, ^s{ui) < implies stability as path perturbations decay in 
time, while $s(uj) < implies an instability as path perturbations are amplified. The 
real and imaginary parts of the complex angular frequency a;, as a function of v : are 



shown in figure 10, We observe that $s(lj) changes sign from positive to negative at 
high velocity v c close to c s . The critical velocity v c is only weakly dependent on the 
Poisson ratio (not shown). Moreover, 3?(cj)^0 at this velocity. Together, these imply a 
high-velocity oscillatory instability with a wavelength A osc that satisfies 

Ksc ~ T d V c ~ t n l . (62) 
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Figure 10. 9(u)) (left) and 5R(a)) (right) as a function of v/c 3 , where Co = ujt^ (see 
|154j for more details). The existence of a linear oscillatory instability is predicted at 
v c ~0.77c s , for which 9(a)) becomes negative with 3?(dj)^0. 

This is an important prediction, suggesting the existence of a high- velocity oscillatory 
instability that is controlled by intrinsic time and length scales associated with the 
breakdown of LEFM near the tip of a crack. The critical velocity v c is about 80% of the 
shear wave-speed, but in light of the various approximations adopted, we do not intend 
to compare its exact value to experimental data. 

6.2. Experimental test 

The theory described in the previous section predicts that, in the presence of a finite 
t n \, causality implies that the singular LEFM fields lag behind the actual tip location 
with a delay time r^oc^. This is linked to a high- velocity oscillatory instability with 
the following properties: 

• The scaled critical velocity for the onset of oscillations v c /c s is close to unity and 
material independent. 
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Figure 11. The oscillation wavelength, A osc changes significantly as the instability 
develops. Shown are typical experiments depicting the wavelength evolution with their 
order of appearance for materials with (a) T = 39 J/m 2 , /i = 168/cPa (b) T = 27 J/m 2 , 
ji = 36kPa. (bottom) Typical photographs of the xy profiles of fracture surfaces 
corresponding to the data sets. 



• The oscillation wavelength X osc is proportional to l n \. 

The best way to test these predictions is to excite the oscillatory instability in 
a variety of different materials, thereby controllably varying both \i and T(v) (and 
hence £ n ^ cf. equation ([53])). These experiments were performed in [99j by preparing a 
variety of polyacrylamide gels, by varying the concentration of the gel components, and 
thereby changing both the elastic moduli and T for each gel composition. In this way we 
constructed gels over a range of shear moduli \i (33</x< 187 kPa) and fracture energies 
24 < T(v c ) < 60 J/m 2 at the critical velocity for the onset of the instability, v c . The 
experiments were performed in thin sheets of dimensions (xxyxz) (130 x 130 x 0.2)mm 
and (200 x 200 x 0.2)mm, where as before x, y and z are, respectively, the propagation, 
loading and thickness directions. These dimensions were large enough so that the 
dynamics at v c occurred in an effectively infinite 2D medium. The sheet thickness was 
small enough (160— 220/iin) to suppress the micro-branching instability to enable single 
crack states to attain velocities to beyond 0.9c s (as in [98]). Once excited, the oscillatory 
wavelength, A osc , may change significantly as the instability evolves, as demonstrated 
in figure [TT] In order to test the theoretically predicted initial wavelength, comparison 
with the theory was performed using only the first excited wavelength. 

In each of the experiments shown, experimental conditions were identical except 
for the material used. The first prediction of the theory is that the scaled critical 
velocity, v c /c S) at the onset of the oscillatory instability is of order unity and material 
independent. This is born out in figure [12^, which demonstrates that v c « 0.9c s for all 
of the materials tested. Figures [12)3 and c demonstrate that variation of the material, 
on the other hand, causes a wide variation of the of the oscillatory wavelength, \ OSC) 
with the elastic modulus of the material. We wish to compare X osc to the non-linear 
scale, £ n i : predicted by the theory. 
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Figure 12. (a) The scaled critical velocity, at the onset of the instability v c ^0.9c s , 
is constant. v c is defined as the maximal velocity prior to the instability onset in 
each material. Symbol colors correspond to the legend of panel b (b) The oscillation 
wavelength, A osc changes significantly with the material, as characterized by fi. (c) 
Typical photographs of the xy profiles of fracture surfaces at the onset of the oscillatory 
instability; from top to bottom: /i= 36, 143, and 168 kPa. Figure adapted from [99] . 



There are a number of ways to estimate or measure t n \. One method is by means 
of equation (53) which states that t n \ oc though with a highly nontrivial pre- 

factor. Using this loading configuration, T{v) could be measured using equation (15), 
by determining Ki(v c ) via the crack-tip curvature (cf. equation (46)), when corrected 
for finite strain [99]. T{y c ) is determined by measuring the crack tip curvature at v c and 
utilizing the universal kinematic functions calculated in the framework of LEFM. These 
functions (see equation 15) become singular as v cr. As a result, the value of T(v c ) 
can be quite sensitive to small uncertainties in cr. 

An additional method for estimating the scale of t n \ is by directly measuring S(v c ) 
in each material. 5(y c ) is a scale that incorporates all of the nonlinear elastic effects in 
the vicinity of the crack tip, and therefore must be proportional to i n \. This method 
has the distinct advantage of using a directly measurable quantity, with relatively small 
uncertainties. 

Comparing the two estimates of i n \ in figure [T3l we see that they are both 



proportional to each other with nearly the same magnitude. The proportionality of 
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Figure 13. A comparison of the two estimates for the nonlinear scale, £ n i, for the 
different materials used (colors correspond to the legend in figure (12)3). We see that 
both estimates, T(v c )//jL and 5(v c ), are linearly dependent with approximately the same 
size. Note that the dashed line, which is a guide to the eye, indicates an offset value 
of ~100-200/xra. 

these two scales, however, is marred by a slight offset of about ~ 100 — 200/im, as 
indicated by the dashed line in the figure. This is perhaps not surprising, since T/fi 
is an estimate of t n \ predicted by the weakly non-linear theory, whereas 5(y c ) is a 
directly measured quantity that accounts for all nonlinear (i.e. deviations from linear 
elasticity) contributions to the crack tip scale (i.e. weak and strong elastic nonlinearities, 
dissipation etc.). 

In figure [14] we compare the wavelength of the first observed oscillation, \ OSC) which 
should correspond to the linearly unstable wavelength in the region of linear growth 
of the instability, to 5(v c ) for the 6 materials used. We indeed find that 5 is directly 
proportional to A osc , as predicted in j!54j . 

In summary, the results of this section conclusively demonstrate that the oscillatory 
instability of fast brittle cracks indeed involves an intrinsic scale that is governed, in 
a large part, by the nonlinear elastic zone surrounding the crack tip. The size of 
this zone quantitatively agrees with the predictions of [154] . These results indicate 
that the nonlinear (and dissipative) zones surrounding the tip of a moving crack are 
not "passive" objects that are simply "dragged along" by the crack tip. Instead, as 
suggested by [74j EHl EQ EH H54j . this region may play an active role in destabilizing 
crack motion. The demonstration of this, as presented in this work, is therefore an 
important step in obtaining a fundamental understanding of the origin of instabilities 
in dynamic fracture. These ideas are as general as the singular behavior that occurs at 
the tip of a moving crack. We believe that it is, therefore, quite likely that dynamics 
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Figure 14. Comparison between the non-linear length-scale S(v = v c ) and the 
oscillation wavelength A osc . Note that the different combinations of /i, T, and e are 
used to produce ~ 15 independent measurements. The dashed line is a guide to the 
eye. Data are taken from [99] . 



of the near-tip zone could play an important role in unraveling the physical mechanism 
driving other instabilities of rapid cracks [59J, EU I10UI 11951 1196] . 

7. Summary and open challenges 

In this article we first briefly reviewed the most well-developed theory of fracture - 
Linear Elastic Fracture Mechanics (LEFM) - and its major predictions. We argued 
that while this theory is very successful in various aspects, it falls short of explaining 
the fast dynamics of a crack once it deviates from a perfectly straight path. Thus, 
high-velocity path instabilities, most notably the side-branching and the oscillatory 
instabilities, remain open problems in this framework. We then summarized recent 
experimental and theoretical progress in understanding the dynamics of rapid brittle 
fracture, with a special focus on instabilities. We have highlighted the central role 
played by near crack front/tip nonlinearities and the associated intrinsic length scales 
in understanding these instabilities. 

From an experimental perspective, we demonstrated that soft brittle elastomers 
mirror in detail the fracture phenomenology of more traditional brittle materials. By 
taking advantage of their significantly reduced wave-speeds we have been able to 
probe the brittle fracture process at length and time scales that were not previously 
directly accessible. These experiments have led to two major outcomes. First, they 



The Dynamics of Rapid Fracture: Instabilities, Nonlinearities and Length Scales 47 

enabled precise quantitative verification of detailed predictions of Linear Elastic Fracture 
Mechanics (LEFM) for the dynamics of straight cracks to an unprecedented degree and 
dynamical range. These experimental tests include both straight cracks propagating in 
an effectively infinite medium and straight cracks propagating in an infinite strip, where 
interactions with external boundaries qualitatively change the nature of the dynamics. 
Second, these experiments have revealed how linear elasticity breaks down near crack 
fronts/tips and clearly indicated the importance of near crack tip nonlinearities and 
the associated intrinsic length scales for crack instabilities, suggesting where essential 
physics is missing. 

From a theoretical perspective, we described in detail the development of a new 
theoretical framework, the weakly nonlinear theory of fracture, which was directly 
motivated by the experimental observations on the breakdown of linear elasticity near 
crack fronts/tips. The basic premise of this theory is that the onset of this breakdown 
process is nonlinear elastic in nature. The predictions of this theory, in particular the 
form of the near crack front/tip singularity and the crack tip opening profile, are in 
excellent quantitative agreement with the experimental measurements. 

This theory also predicts the existence of an intrinsic, i.e. geometry and loading 
independent, length scale that emerges from a competition between linear elastic and 
weakly nonlinear elastic deformations near crack fronts/tips. Based on the existence of 
this finite length scale, an equation of motion for the direction of crack propagation in 
2D has been described. This extension of the nonlinear theory to path dynamics predicts 
a high-speed oscillatory instability whose wave-length is determined the intrinsic length 
scale. We have shown that this prediction is supported by experiments on a variety of 
different soft brittle materials. In our view, this is an important result that demonstrates 
that the near crack front /tip region not only accounts for the dissipation that accompany 
crack propagation, but may also play a central role in determining crack stability. 

The work described in this paper is a basic attempt to account for elastic nonlinear 
response of materials in the vicinity of a crack tip. We have shown that new physical 
effects result when only weak elastic nonlinearity is taken into account. What happens 
when a material is so tough as to enable strong elastic nonlinearities to occur on large 
scales? One such example is the fracture of rubber, where strains larger than unity are 
easily obtained. When such enormous energy densities are reached, it is possible that an 
underlying premise of fracture mechanics, that energy must be transported from remote 
distances to enable fracture, should be re-examined. Work in this direction has, for 
example, demonstrated that in such circumstances supersonic tensile fracture can take 
place [T02irT26irT97] . 

The oscillatory instability discussed in detail in this review has been observed in a 
range of soft brittle materials, where cracks have been driven to extremely high velocities 
and in which elastic nonlinearities are both pronounced and spatially well-separated from 
dissipative nonlinearities. It remains to be seen whether it can be observed in more 
"standard" brittle materials, where no clear separation exists between nonlinear elastic 
scales and dissipative ones, in spite of the enormous technical difficulties involved. One 
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may speculate that the existence of t\ n , whether associated with elastic or dissipative 
nonlinearities, is sufficient for the existence of the oscillatory instability. This challenge 
might be at least partially addressed through advanced numerical simulations. 

Within a broader context, we believe that an additional crucial step in pushing the 
field forward is the substantiation and extension of the 2D equation of crack propagation 
and the development of the 3D counterpart for crack fronts. Addressing this challenge 
from both theoretical and experimental perspectives, and making progress along these 
lines, may open the way for more systematic path stability analysis and enhance our 
ability to quantitatively predict the failure dynamics of materials and interfaces. This 
may have important consequences for various problems in a wide range of fields, ranging 
from materials science to biology. 

Such progress may also pave the way to resolve one of the most resistant puzzles 
in the field of dynamic fracture - the side-branching instability. As was shown above, 
aspects of this instability possess an intrinsic 3D nature. One might then surmise 
that to "crack" this instability one might have to further extend our understanding 
of 3D crack front dynamics pPfl 11471 I198j . The latter is intimately related to other 
open problems in the field, such as the roughness of crack surfaces JT91 EQl 11991 [200J, 
crack front waves [T71 HHJ 12011 1202j and the stability of crack fronts to an ensemble of 
continuous perturbations [22l 12031 1204j . 
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